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^ : ABSTRACT 



We present a 5-8 /im analysis of the Spiteer-IRS spectra of 71 ultraluminous 



■ infrared galaxies (ULIRGs) with redshift z < 0.15, devoted to the study of 
O ' 

I the role of active galactic nuclei (AGN) and starbursts (SB) as the power 

source of the extreme infrared emission. Around ^5 /im an AGN is much 
brighter (by a factor w30) than a starburst of equal bolometric luminosity. 

X ■ 

■ This allows us to detect the presence of even faint accretion-driven cores 

inside ULIRGs: signatures of AGN activity are found in ~70 per cent of our 
sample (50/71 sources). Through a simple analytical model we are also able 
to obtain a quantitative estimate of the AGN/SB contribution to the overall 
energy output of each source. Although the main fraction of ULIRG luminosity 
is confirmed to arise from star formation events, the AGN contribution is 
non-negligible ('-~^23 per cent) and is shown to increase with luminosity. The 
existence of a rather heterogeneous pattern in the composition and geometrical 
structure of the dust among ULIRGs is newly supported by the comparison 
between individual absorption features and continuum extinction. 
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1 INTRODUCTION 



Since a significant population of ultraluminous infrared galaxies (ULIRGs, i.e. the sources 
whose 8-1000 /xm luminosity exceeds > IO^^Lq) was discovered during the IRAS all- 
sky survey (Houck et al. 1985), the origin of their huge infrared (IR) emission has been 
widely debated. Extreme IR activity is known to be closely related to interacting or merging 
systems: optical and near-IR imaging shows that a large part of ULIRGs is indeed result- 
ing from the collision of gas- rich galaxies (e.g. Clements et al. 1996; Scoville et al. 2000). 
Multiwavelength observations have pointed out that the typical spectral energy distribution 
(SED) of ULIRGs is characterized by a prominent far-IR peak due to dust reprocessing of 
higher-frequency primary radiation, whose elusive source is supposed to be a combination 
of merger-induced bursts of star formation and gas accretion onto a supermassive black 
hole. Despite their low frequency and limited contribution to the energy density of the local 
Universe (Soifer & Neugebauer 1991), nearby ULIRGs represent the counterparts of both 
the submillimetre sources (Blain et al. 2002, and references therein) and the 24 //m-selected 
galaxies that are responsible for the dominant energy output at high redshift (Caputi et 
al. 2007). Understanding their nature and physical properties is therefore a task of great 
importance, with far-reaching implications for models of structure formation and galaxy 
evolution. In particular, many efforts have been made in order to assess the contribution 
of the starburst (SB) and active galactic nucleus (AGN) components to ULIRGs: such an 
estimate would give deeper insight into the history of star formation over cosmic time and 
the origin of the X-ray and far-IR backgrounds. 

All the diagnostic studies concerning ULIRGs have to come to terms with the great opacity 
of their nuclear regions, that precludes a straight identification of the hidden power source. 
Nevertheless the mid-IR thermal emission is expected to suffer from little extinction, and at 
the same time to preserve some relic signatures of the underlying engine, making this wave- 
length range very attractive for solid ULIRG diagnostics. With the coming of the Spitzer 
Space Telescope (Werner et al. 2004) and its Infrared Spectrograph (IRS; Houck et al. 2004), 
high-quality mid-IR spectra of several tens of local ULIRGs have become available. We have 
therefore been able to perform a 5-8 /xm spectral analysis of a sample of 71 local ULIRGs, 
separating the AGN and SB components in the observed emission and providing a quantita- 
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tive estimate of the AGN/SB contribution to the bolometric luminosity of each sourcqj. The 
prehminary results have been anticipated in a recent letter (Nardini et al. 2008, hereafter 
Paper I), mainly focused on the presentation of the decomposition method. In the present 
work we will discuss in detail our findings and investigate their implications. The outline of 
this paper is the following: in Section [2] we briefly review the main studies, from radio to 
X-rays, that have contributed to clarify the multiwavelength picture of ULIRGs. Section [3] 
illustrates the physical motivations for 5-8 fim diagnostics and describes how recent Spitzer 
observations have prompted our AGN/SB decomposition effort. The ULIRG sample is pre- 
sented in Section 4, while Section 5 concerns observations and data reduction. In Section 6 
we summarize the fitting model and all the analytical steps leading to the estimate of the 
AGN/SB contribution to the bolometric luminosity. The comparison between our results and 
those obtained at different wavelengths is provided in Section 7, and the main implications 
of our findings are discussed in Section 8. The conclusions are drawn up in Section 9, and 
all the fits are shown in the Appendix along with further considerations on some notable 
objects. 



2 REVIEW OF ULIRG DIAGNOSTICS 

The detection of a faint or obscured AGN and the estimate of its contribution to the bolo- 
metric luminosity of the hosting galaxy are the main open issues in the study of ULIRGs. 
This task is complicated because of the dust reprocessing, that smooths the intrinsic dif- 
ferences between the primary AGN and SB radiation fields. A detailed knowledge of the 
broadband intrinsic emission of both active galactic nuclei and starbursts is then essential 
in order to find some signatures of the buried energy source. 

In the optical ULIRGs are classified on the basis of intensity ratios, involving emission lines 
with different excitation energy to test the hardness of the radiation field. However, dust 
extinction and reddening affect the use and reliability of such optical diagnostics in various 
ULIRGs; moreover the possibility of a differential extinction of the AGN and SB components 
is difficult to take into account. Qualitatively, compelling evidence of AGN activity is found 
in a modest fraction of the 118 ULIRGs included in the IRAS 1 Jy sample (Kim & Sanders 
1998), with a possible trend of increasing AGN detection rate in sources of higher luminosity 



The bolometric luminosity of a ULIRG is usually dominated by its 8—1000 fim emission: throughout this paper the terms 
bolometric luminosity and IR luminosity are consequently used with an equivalent meaning. 
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(Veilleux, Kim & Sanders 1999). At Lm > 2 x lO^^L© nearly a half of IR galaxies exhibit the 
typical features of Seyfert galaxies, and the fraction of AGN-like objects is even larger if one 
includes the sources classified as Low-Ionization Nuclear Emission-line Regions (LINERs). 
The connection between AGN and LINERs is anyway controversial, because the ionization 
degree of the latter can be a possible consequence of the supernova-driven shocks and the 
violent galactic winds pervading a starburst as well (Taniguchi et al. 1999; Lutz, Veilleux & 
Genzel 1999). This leaves unsolved the issue about the fraction of ULIRGs actually hosting 
an AGN. 

The dust opacity affecting the optical, ultraviolet (UV) and near-IR domains turns irrele- 
vant at radio wavelengths, which in principle also benefit of a much higher angular resolution 
accessible through interferometry. Unfortunately these advantages are reduced by the low 
radio emission of ULIRGs. Nevertheless, although jet-hke structures or extended lobes have 
never been found in ULIRGs, many of them present compact radio cores whose power does 
not correlate with the Ha luminosity nor with other usual indicators of star formation (Na- 
gar et al. 2003). Better spatial resolution is required to confirm the AGN nature of such 
sources and rule out any hypothesis of a compact cluster of stars. 

Speculatively, the optimal spectral region to clear all the ambiguities and detect the AGN 
component is the 2-10 keV band, where it is possible to look through the obscuring material 
and take advantage of the sharp difference between the hard X-ray features and luminosity 
of starbursts and AGN. As a matter of fact, the circumnuclear mixture of gas and dust 
can be Compton-thick, with equivalent column densities in excess of 10^^ cm~^, and bury an 
AGN even in this energy range. Moreover ULIRGs are faint emitters, since their typical hard 
X-ray flux is usually far less than ~10~^ times their IR flux (e.g. Risaliti et al. 2000). Hence 
the current instrumentation allows the X-ray diagnostics for a handful of bright sources only: 
intense star formation activity is confirmed to be the dominant power supply for extreme IR 
luminosity, and even if convincing traces of elusive AGN (i.e. missed in the optical, Maiolino 
et al. 2003) are found, the conclusions are still tentative (Franceschini et al. 2003; Ptak et 
al. 2003). 

In the mid- and far-IR, the reprocessed radiation still provides some clues to the enigmatic 
energy source, that can be recognized e.g. in the circumnuclear dust composition, in the SED 

large-scale profile and/or in the appearance of coronal lines. An empirical indicator of star 
formation is usually considered to be the group of emission features centred at 3.3, 6.2, 7.7, 
8.6 and 11.3 /im, and attributed to polycyclic aromatic hydrocarbons (PAHs): these features 
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are prominent in dusty environments exposed to soft UV radiation but strongly suppressed 
in the presence of harder radiation fields, probably because of the alteration/destruction of 
their carriers (Voit 1992). The concurrent use of PAH strength and appropriate line ratios, 
e.g. [O lv]/[Ne ll] and [Ne v]/[Ne ll], results in an effective probe into the nature of ULIRGs 
(Genzel et al. 1998). Yet the diagnostic methods based exclusively on emission lines can 
not be easily applied to faint sources. This involves the necessity of introducing the contin- 
uum properties and mid-IR colours as alternative diagnostic tools, as suggested by Laurent 
et al. (2000), with the replacement of the line ratio with the ratio between the warm and 
the hot dust continua. A critical shortcoming actually remains, consisting in the possible 
non-detection of a weak AGN inside the ULIRG: an unresolved AGN component which is 
heavily absorbed or intrinsically feeble can vanish within the much more diffuse and intense 
starburst emission. 

3 THE AGN/SB SPECTRAL DECOMPOSITION AT 5-8 

In order to overcome the potential AGN-dilution effect, it is necessary to select the wave- 
length range in which the imbalance between the AGN and SB components is maximum. 
By examining the average SEDs of both AGN and starbursts one immediately evinces that 
an AGN is much brighter than a starburst with equal bolometric luminosity not only in the 
hard X-ray domain, but also in a narrow IR band stretching from ~3 fim, where the red tail 
of direct stellar emission begins to vanish, up to ~8 /im, at the onset of the possible silicate 
absorption feature centred at 9.7 /im (see Fig{T]). This difference can be ascribed to the more 
intense UV radiation field of the AGN, which is able to heat up to the sublimation temper- 
ature the dust located in proximity of the central energy source {Tgub ~1500 K for graphite 
grains). Incidentally, both the above-mentioned mid-IR colour diagnostic and the previous 
identification of a subset of warm ULIRGs (Sanders et al. 1988) are based on a completely 
similar argument. It is only in the recent years that this spectral region and its potential in 
the study of ULIRGs have been exploited quantitatively (e.g. Risaliti et al. 2003, 2006a). In 
particular, L-band spectroscopy (~3-4 fim) has proven to be a powerful method to uncover 
a buried AGN inside a ULIRG, thanks to diagnostic features like the equivalent width of the 
3.3 /im PAH emission and the continuum slope (Risaliti et al. 2006b). Specifically, a steeply 
rising continuum hints at the reddening of a point-like energy source; an obscured AGN 
interpretation can be also supported by the simultaneous detection of deep absorption pro- 
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Figure 1. Dependence from wavelength in tlie Spitzer-IRS operative window of the AGN over SB brightness ratio for 
equal bolometric energy output. The curves are obtained by dividing the observed emission of IRAS 07598+6508 (red line), 
IRAS 08572+3915 {blue line) and IRAS 12127—1412 [green line) by that of our SB template, whose definition is described 
later on in this Section. The shaded regions show the la rms uncertainty, and the spectral resolution has been degraded in 
order to smooth the dependence of the ratios from narrow features. The three objects appear as AGN-dominated at a visual 
inspection, with moderate reddening, significant absorption and/or suppression of aromatic features. After our analysis, the 
AGN bolometric contribution to each source turns out to exceed 80 per cent. It is clear how the AGN spectral enhancement 
rapidly decreases with increasing wavelengths, giving rise to a dilution effect that precludes the detection of intrinsically faint 
or heavily obscured AGN components. Concerning the red curve, the prominent 10 ^m bump and the following plateau are 
exclusively due to the lack of silicate absorption, which is unusual in composite ULIRGs. 



files at ~3.4 and ~4.5 /im due to aliphatic hydrocarbons and carbon monoxide, respectively 
(Sani et al. 2008). Ground-based observations are anyway affected by the strong variability 
of atmospheric transmission and the high thermal background. For this reason the quality 
of the current L- and M-band spectra, even if obtained with large-aperture facilities like 
SUBARU and VLT, allows the determination of the AGN and SB components for only the 
~20 brightest ULIRGs. In this view, the spectral capabilities at 2.5-5 yum of the AKARI 
satellite look very promising (e.g. Oyabu et al. 2009). 

A significant progress has been possible thanks to Spitzer, whose spectrograph operates long- 
ward of 5 /im and provides the natural extension of our diagnostic study. The luminosity 
enhancement of the AGN over the SB component is expected to substantially fade and turn 
over with increasing wavelength (FigJT]), none the less at 5-8 fim the considerably higher 
accuracy achieved from space still overrides this effect and allows a precise determination 
of the AGN and starburst properties in a composite ULIRG. Before proceeding with the 
details of our diagnostic technique, it is worth exploring the possibility for a luminous but 
extremely obscured AGN to be missed even at 5-8 fim. We consider the test case of a source 
whose IR luminosity arises in equal parts from nuclear activity and star formation. A con- 
tribution of 20 per cent to the observed 5-8 fim emission can be assumed as a reasonable 
detection limit for the AGN component, which intrinsically is expected to outshine the star- 
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burst counterpart by a factor ~20-30 in this spectral range. Consequently, the AGN should 
undergo a flux attenuation of ~2 orders of magnitude in order to fall below the detection 
threshold. In terms of a standard extinction law this corresponds to a screen extinction of 
Ay ~200 mag {AQ^rn/Ay ~0.025; Nishiyama et al. 2008, 2009). This scenario is possible but 
unlikely (see also Section [6] and the discussion in Genzel et al. 1998). We are then confident 
that virtually all the luminous AGN components inside ULIRGs can be probed through 
5-8 yum Spitzer-YR^ spectra. 

The unprecedented sensitivity of the mid-IR capabilities onboard Spitzer, in fact, has made 
it possible to reconstruct with high precision the SEDs of large samples of luminous active 
nuclei and starburst galaxies, and a surprising spectral homogeneity has been found within 
the two separate classes. At the highest luminosities this is especially true in the 5-8 fim 
range, with nearly constant shapes of both the aromatic complex in starbursts (Brandl et al. 
2006, hereafter B06) and the continuum in AGN (Netzer et al. 2007, hereafter N07). Larger 
differences come out at longer wavelengths, where the temperature and emissivity of the 
different dust components as well as radiative transfer need to be taken into account. Our 
strategy then aims at disentangling the AGN and SB contributions to the observed 5-8 fim 
emission of ULIRGs by means of spectral templates, whose success hinges upon the large 
difference between the average AGN/SB emission and the small dispersion within the AGN 
and SB spectral classes. In the following we explain in detail the physical reasons which 
justify the selection of the AGN and SB templates for this work. 

Starburst. As previously mentioned, the mid-IR spectra of starburst galaxies display lit- 
tle variations from one object to another, especially in the 5-8 /im wavelength range. At first 
sight this appears to be a singular circumstance, because of the large number of physical 
variables involved in determining the observed properties of a starburst: these include the 
initial mass function, the duration and evolutionary stage of the single bursts, the properties 
of the dust grains and those of the interstellar environment, and obviously the total mass 
and luminosity of the system. It has already been suggested in B06 that such a remark- 
able similarity can derive from the spatial integration over a huge number of unresolved 
star-forming spots: this is indeed an interesting point that should be investigated further. 
The typical SB emission in the 5-8 fim interval is characterized by two prominent aromatic 
features centred at 6.2 and 7.7 fim. In SB-dominated ULIRGs the shape of these lines and 
their intensity ratio slightly differ from those observed in lower luminosity starbursts (like 
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NGC 7714), probably because of higher obscuration in the former (Rigopoulou et al. 1999). 
Moreover the reliability of PAHs as quantitative tracers of star formation is still debated 
(Peeters, Spoon & Tielens 2004; Forster Schreiber et al. 2004). In order to properly de- 
scribe the SB component in our spectra, and check if the prescription of little dispersion 
may be applied in the ULIRG luminosity range as well, we performed a preliminary analysis 
on a representative subset of sources that other multiband diagnostics have established to 
host powerful starbursts. The 5-8 fim emission of these sources has been approximated by 
means of two gaussian profiles for the aromatic features and a power-law continuu As a 
result, the main properties of the PAH features turn out to be the same observed in lower 
luminosity systems. Their equivalent width is confirmed not to depend on the global IR 
output, and the modest spread concerning their relative strength has already been related 
to internal obscuration; the slope of the underlying continuum is F^f, >4, roughly uniform 
and well above the value assumed for the AGN component (see later). This tight correlation 
between the properties of the PAH features and the dust continuum is brought out also in 
B06, hence we are justified in the use of a spectral template to represent through a single 
scale factor the SB component of ULIRGs. A complete characterization of the whole SB 
spectrum would of course require a detailed model involving many parameters, and this 
indeed remains unavoidable when considering larger spectral domains. 

Instead of adopting a prototypal SB template from the literature, we have constructed a 
new one as the average spectrum of the five brightest SB-dominated sources in our sam- 
ple, i.e. IRAS 10190+1322, IRAS 12112+0305, IRAS 17208-0014, IRAS 20414-1651 and 
IRAS 22491 — 1808. The resulting slope is F^fc ~4.2, and the equivalent widths of the aro- 
matic features are 0.43 and 0.88 ^m: such values do not actually differ from the measures 
of B06, as can be easily evinced from Figj2l The comparison between our SB template and 
that of B06 over the entire IRS spectral range has already been shown in Paper I, in order 
to emphasize their divergence just longward of ~8 fim. This is anyway displayed again in 
the box of FiglSl where we focus instead on their similarity over the wavelengths of inter- 
est. Such similarity (and the fact that the PAH equivalent widths are the largest among 
ULIRGs) rules out any possibility of significant AGN contamination in our SB template- 



^ Throughout this paper we describe the emission in terms of the flux density fv(^), hence all the spectral indices of the 
power-law components in the following are defined as /j/(A) oc A'". 

^ We also note that in this case our AGN/SB spectral decomposition would have brought out the necessity of a negative AGN 
contribution in order to optimally fit the shape of a genuine SB-dominated ULIRG. 
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Figure 2. Comparison between our SB template (blue solid line) and tiiat of B06 {deep green dashed line) at 5-8 fim. Their 
beiiaviour up to ~30 fim is outlined in a log-log scale in the top left corner. The single spectra of the five sources employed 
in our template have also been plotted with lines of different colour: IRAS 10190+1322 {light green), IRAS 12112+0305 {red), 
IRAS 17208-0014 {magenta), IRAS 20414-1651 {orange) and IRAS 22491-1808 {cyan). 

Summarizing, the spectral slope of the observed 5-8 fim ULIRG continuum is very sensitive 
to the possible AGN contribution. The PAH features and their intensities play of course an 
important role in modelling the spectral shape, but the dispersion in their properties which 
is found in pure starbursts is actually too small to alter the determination of the continuum, 
mimicking the presence of a bolometrically significant AGN in a composite source. 

Active galactic nucleus. At the wavelengths of concern the AGN emission is dominated 
by the cooling of small dust grains that are transiently heated up to temperatures close to 
the sublimation limit. Such a process is expected to produce a nearly featureless power-law 
continuum, and yet aromatic features are found in several quasars, as confirmed by recent 
observations (Schweitzer et al. 2006; Lutz et al. 2007, 2008). The inferred star formation rate 
can be very high, and raises the question on the extent of the stellar contamination in the 
observed spectra of quasars and the reliability of a pure AGN template. In N07 the far-IR 
emission of a large set of Palomar-Green (PG) quasars is ascribed to cold dust in extended 
regions, and then suspected of having a stellar origin; the contribution of star formation to 
be subtracted in order to determine the intrinsic SED of a type 1 AGN is obtained from the 
properties of 12 SB-dominated ULIRGs. The resulting average spectrum can be modelled 
by a single power law from ~3 fim to the 9.7 fim silicate bump (which usually appears 
in emission whenever the source is not obscured). The most striking outcome is that the 
spectral dispersion at these wavelengths is very small if compared to the large differences 
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Figure 3. The blue solid line represents the average L-band spectrum of a source whose near-IR emission is dominated by an 
unobscured AGN; the shaded area encompasses the 1(T rms dispersion. The best fit {red dashed curve) requires two components; 
a power law with spectral index i.5f±0.09, and a gaussian line with equivalent width < 3 nm, far below the typical values of 
the PAH emission in starburst galaxies (~100 nm). 

found in the far-IR. It is then reasonable to adopt a power-law component to reproduce the 
5-8 fim AGN flux density. The choice of the spectral index however is not immediate, and 
has to be briefly discussed. A coarse estimate, as deduced from the SEDs provided by N07, 
yields Tagn ~0.7; a better value can be anyway achieved from our L-band study. The reason 
is that we are dealing with an active nucleus component in a ULIRG environment, and 
the overall dust properties are expected to differ from those of the PG quasars in the N07 
sample. In FigJSlwe show the spectral template obtained by averaging the L-band emission 
of four ULIRGs in our sample known to host a fairly unobscured active nucleus: MRK 231, 
IRAS 15462+0450 and IRAS 21219-1757, which are optically classified as type 1 Seyfert 
galaxies, and IRAS 05189—2524, which is a type 2. For all these sources the L-band emission 
consists of a strong continuum arising from the dust heated by the AGN, as also confirmed 
by the small equivalent width of the 3.3 /im aromatic feature, only reaching a few nm while 
usually exceeding ~100 nm in pure starbursts (Moorwood 1986; Imanishi & Dudley 2000). 
The mean value of the continuum slope is ~1.5, and such a value has then been assumed 
for the present analysis. Its extrapolation to 5-8 fim is supported by N07, since no signifi- 
cant change in the continuum gradient is expected. Moreover the same 5-8 fim spectra of 
the four template sources are fully consistent with the stretching of their L-band slope, in 
spite of the growing aromatic emission. The slight discrepancy between these estimates of 
the AGN spectral index, based on different samples, seems to be a natural consequence of 
the greater obscuration characterizing the ULIRG environment, rather than of a significant 
stellar contamination even at such short wavelengths. As mentioned, the unobscured AGN 
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components of this work and the PG quasars of N07 (even though intense star formation 
is in place) are different populations, in terms of both the luminosity of the AGN and the 
properties of the absorbing medium. A higher covering factor of the dust surrounding the 
AGN component in a ULIRG environment can account for a steeper intrinsic slope. The 
covering factor is indeed shown to decrease as a function of the AGN luminosity (Maiolino 
et al. 2007). 

Beyond these aspects, another crucial effect has to be considered: the active nucleus compo- 
nent is still extremely compact, and the AGN-related hot dust emission can be extinguished 
by an external screen, consisting either of the outer layers of the obscuring torus predicted in 
unification models (Antonucci 1993) or of a thick absorbing cloud along the line of sight. We 
have therefore introduced a wavelength dependent attenuation factor e"'^^'^-', assuming for 
the optical depth a power-law behaviour r(A) cx A^^''^^ (Draine 1989). On the contrary, the 
SB component can easily extend over several hundreds of pc and therefore can be affected 
exclusively by internal extinction, which is already embodied in the shape of the template: 
in fact, even if the SB extinction becomes apparent when considering the silicate absorption 
profiles at longer wavelengths, the 5-8 /xm emission does not seem to be significantly modi- 
fied except for the limited dispersion in the PAH ratio. 

In Section [8] we will deal again (from a quantitative point of view) with the intrinsic slope 
of the AGN continuum emission, investigating its interplay with the shape of the extinction 
law, and how any change can affect the AGN/SB spectral decomposition and the subsequent 
estimate of the relative AGN/SB bolometric contribution. 

4 THE ULIRG SAMPLE 

The present sample consists of 71 objects and is fully representative of the ULIRG population 
in the local Universe. In order to avoid a bias towards AGN activity, the flux limitation 
has been adopted at 60 fim. The 60 fim flux density is in fact a good proxy of the cold 
dust component, and this translates into a fairly unbiased selection with respect to the 
nature of the energy source. All the ULIRGs in the already mentioned 1 Jy sample with 
available Spitzer-IKS observations and the additional requirement of a redshift z < 0.15 
have been selected, for a total of 63 objects. With respect to Paper I, three extra sources 
have been added, i.e. IRAS 02021-2103, IRAS 07598+6508 and IRAS 08559+1053. The 
broad absorption line (BAL) Seyfert 1 galaxy IRAS 07598+6508, that was initially dropped 
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Table 1. General properties of the 71 sources in our sample, z, D^: Redshift and luminosity distance (in Mpc). /12, /2B, 
/so, /lOo: IRAS flux densities (in Jy, with the colon indicating the revisions in Sanders et al. 2003). L/zj: Total IR luminosity 
computed according to Sanders & Mirabel (1996). Class: Morphological classification as defined in Veilleux, Kim & Sanders 
(2002). Ilia: Wide binary (apparent separation > 10 kpc); Illb: Close binary (apparent separation < 10 kpc); IVa: Diffuse 
merger; IVb: Compact merger; V: Old merger. The entries within brackets are derived from Due, Mirabel & Maza (1997) and 
Scoville et al. (2000). 
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since a significant fraction of its mid-IR emission was expected to liave a non-thermal origin, 
and IRAS 08559+1053, tliat was only observed during the in-flight calibration phase of 
Spitzer (by J.R.Houck), have been retrieved for completeness. We point out that the volume 
restriction is absolutely arbitrary in this context, since it is directly imported from our 
previous studies, being it necessary to apply our L-band diagnostics. The 1 Jy sample indeed 
contains 70 objects meeting this selection requirement. No observation is anyway available 
for the seven missing sources: only four has been recently observed with Spitzer, but archive 
data are not provided as yet; the other three have not been observed. Within this uncovered 
subset five objects are optically classified as Seyfert 2 galaxies, while two are unclassified. The 
consequent bias is then minor, and will be discussed later on when we assess the incidence 
of black hole accretion and star formation at extreme IR luminosities. 
Summarizing, all the selected targets obey the following criteria: (1) a. redshift z < 0.15 
and (2) a 60 fim IRAS flux density /go > 1 Jy. The statistics is already large enough 
to derive general conclusions about the local ULIRG population. However, being the 1 Jy 
sample also characterized by constraints on the position of the target, a systematic search 
for ULIRGs inside the IRAS 2 Jy or 1.2 Jy all-sky samples (Strauss et al. 1992; Fisher et al. 
1995, respectively) would yield a wealth of additional sources fulfllling the present criteria. 
We decided to complete our sample by adding only the eight sources drawn from the 2 Jy 
sample which are already present in Genzel et al. (1998), since they are among the brightest 
and best-studied ULIRGs. 

The main properties of all the sources in our sample are listed in TabUl NGC 6240 is the only 
source that does not strictly comply with the ULIRG classiflcation, because of its slightly 
lower IR luminosity. It is anyway an advanced merger showing all the morphological and 
physical properties of ULIRGs, and is usually included in this class. Concerning this, it is 
also important to specify that the total IR luminosity of our sources has been computed 
according to the broadband flux equation in Sanders & Mirabel (1996): 



where Fj^ is the total IR flux in units of ergs s~^ cm~^ and fu, /25, /eo and /loo are the 
IRAS flux densities in Jy. When only an upper limit is available for /12 (and 725) we have 
simply assumed the one half of this limit. Such strategy provides a good approximation of 
the true flux density if compared to the few estimates proposed in Kim & Sanders (1998). 



Fjn = 1.8 X 10-1^(13.48/12 + 5.I6/25 + 2.58/60 + /loo), 



(1) 
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Table 2. Observation log of our ULIRG sample. All the sources except one have been observed within the following programs: 
Spectroscopic study of distant ULIRGs II (PID 105, PI J.R.Houck), Buried AGN in ultraluminous infrared galaxies (PID 2306, 
PI M.Imanishi), The evolution of activity in massive gas-rich mergers (PID 3187, PI S.Veilleux). TsL- Integration time in 
number of cycles times seconds for the IRS Short-Low orders (SL2: ~5.2-7.7 /im; SLl: ~7.4— 14.5 fim). Since each cycle consists 
of two nod positions, the total observing time per slit is obtained by multiplying these entries by an additional factor 2. 



Object 


PID 


Date (UT) 


Tst 


Object 


PID 


Date (UT) 




ARP 220 


105 




3x14 




2306 


9nn/i Till 1 

ZUWi JUl ID 


2x60 


IRAS 


00091-0738 


3187 


ZUUO J nil ou 


2 X 60 


IRAQ 1zL1Q7^nS1^ 


3187 


9nn^ Tn'eiK 1 ^ 
ZUUO Cclj LO 


2 x 60 


IRAS 


00188-0856 


105 


2003 Dec 17 


2 X 60 


IRAS 14252 — 1550 


2306 


2004 Jul 17 


2 x 60 


IRAS 


00456-2904 


3187 


2005 Jul 14 


2 X 60 


IRAS 14348 — 1447 


105 


2004 Feb 07 


1 x 60 


IRAS 


00482-2721 


3187 


2005 Jul 07 


2 X 60 


IRAS 15130—1958 


3187 


2005 Mar 15 


2 X 60 


IRAS 


01003-2238 


105 


2004 Jan 04 


1 X 60 


IRAS 15206+3342 


105 


2004 Jun 24 


1 x 60 


IRAS 


01166-0844" 


3187 


2005 Jan 03 


2 X 60 


IRAS 15225+2350 


2306 


2005 Feb 07 


2 X 60 


IRAS 


01298-0744 


105 


2005 Jul 14 


2 X 60 


IRAS 1 595n^^finQ 


105 


2004 Mar 04 


3x14 


IRAS 


01569-2939 


2306 


2004 Jul 18 


2 X 60 


IRAS 15462—0450 


105 


2004 Mar 02 


1 x 60 


IRAS 


02021-2103 


3187 


2005 Jan 15 


2 X 60 


TR A S 16090 — 0139 


105 


9nn^ Alio- 


1 x 60 


IRAS 


02411+0353 


2306 


9nn^ Tan ^ A 

iUUO J an L'i 


2 X 60 


i Xv^X O i- D i- O D ^ U -L ^ D 


3187 


ZUUO ivi ai LO 


2 X 60 


IRAS 


03250+1606 


3187 




2 X 60 




2306 


9nnA Till 1 A 

ZUU^ J Hi -L^ 


4x 60 


IRAS 


04103-2838 


3187 




2 X 60 




2306 


onnA Till lA 

ZUU^ JUl ±1 


2 X 60 


IRAS 


05189-2524 


105 


9nn4 A/Tar- 99 
zuu^ iviai zz 


3x14 


TR A Q 1 f^ARVJ-'^AA'? 


2306 


9nn4 Till 1 7 

ZUU^ JUl LI 


2 X 60 


IRAS 


07598+6508 


105 


2004 Feb 29 


3x14 


IRAS 17028+5817 


2306 


2004 Jul 17 


2x60 


IRAS 


08559+1053' 


666 


2003 Nov 23 


2x60 


IRAS 17044+6720 


2306 


2004 Jul 17 


2x60 


IRAS 


08572+3915 


105 


2004 Apr 15 


3x14 


IRAS 17179+5444 


105 


2004 Apr 17 


2x60 


IRAS 


09039+0503 


3187 


2005 Apr 18 


2x60 


IRAS 17208-0014 


105 


2004 Mar 27 


3x14 


IRAS 


09116+0334 


2306 


2005 Apr 21 


2x60 


IRAS 19254-7245 


105 


2005 May 30 


3x14 


IRAS 


09539+0857 


3187 


2005 Jun 05 


2x60 


IRAS 20100-4156 


105 


2004 Apr 13 


1x60 


IRAS 


10190+1322'^ 


3187 


2005 May 22 


4x60 


IRAS 20414-1651 


105 


2004 May 14 


1x60 


IRAS 


10378+1109 


105 


2005 Jun 08 


2x60 


IRAS 20551-4250 


105 


2004 May 14 


2x14 


IRAS 


10485-1447 


3187 


2005 May 23 


2x60 


IRAS 21208-0519^* 


3187 


2004 Nov 13 


2x60 


IRAS 


10494+4424 


2306 


2004 Nov 17 


2x60 


IRAS 21219-1757 


3187 


2004 Nov 16 


2x60 


IRAS 


11095-0238 


105 


2005 Jun 07 


2x60 


IRAS 21329-2346 


3187 


2004 Nov 16 


2x60 


IRAS 


11130-2659 


2306 


2005 Jul 12 


2x60 


IRAS 22206-2715 


3187 


2004 Nov 15 


2x60 


IRAS 


11387+4116 


2306 


2005 Jan 11 


2x60 


IRAS 22491-1808 


105 


2004 Jun 24 


1x60 


IRAS 


11506+1331 


3187 


2005 May 25 


2x60 


IRAS 23128-5919 


105 


2004 May 11 


3x14 


IRAS 


12072-0444 


105 


2004 Jan 06 


1x60 


IRAS 23234+0946 


3187 


2004 Dec 13 


2x60 


IRAS 


12112+0305 


105 


2004 Jan 04 


3x14 


IRAS 23327+2913 


2306 


2004 Dec 08 


2x60 


IRAS 


12127-1412 


3187 


2005 Jun 30 


2x60 


MRK 231 


105 


2004 Apr 14 


2x14 


IRAS 


12359-0725 


2306 


2005 Jun 30 


2x60 


MRK 273 


105 


2004 Apr 14 


2x14 


IRAS 


13335-2612 


3187 


2005 Feb 15 


2x60 


NGC 6240 


105 


2004 Mar 04 


2x14 


IRAS 


13454-2956^" 


3187 


2005 Jul 14 


4x60 


40 +12.50 


105 


2004 Jan 07 


3x14 


IRAS 


13509+0442 


2306 


2004 Jul 17 


2x60 


UGC 5101 


105 


2004 Mar 23 


3x14 


IRAS 


13539+2920 


2306 


2005 Feb 07 


2x60 











Notes. " Southern nucleus. The northern one is also observed but is much fainter. ^ Observation performed during the 
in-orbit checkout phase (IRS campaign P: PID 666, PI J.R.Houck). " Both nuclei are observed but apparently not completely 
resolved. Northern nucleus. The southern one is also observed but is much fainter. 

The luminosity distances have been obtained making use of a standard cosmology with 
Ho = 70 km s-i Mpc-\ fi^ = 0.27 and = 0.73 (Hinshaw et al. 2009). 



5 OBSERVATIONS AND DATA REDUCTION 

The spectroscopic observations of all the sources in our sample were obtained with the 
Spitzer-IRS low-resolution modules, within three different programs devoted to buried ac- 
tivity in ULIRGs and massive merger^l: the program IDs are 105, 2306 and 3187, and the 
principal investigators are J.R.Houck, M.Imanishi and S.Veilleux, respectively. The details 
concerning each observation are listed in Tabj2l The local ULIRGs of our sample are bright 



* The only exception is IRAS 08559+1053, as mentioned. 
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objects, and their signal can be easily distinguished from the background emission, mostly 
due to zodiacal light. After having verified that a refined reduction of the data was not 
necessary for our purposes, we decided to perform our analysis beginning from the coadded 
images provided by the Spitzer Science Center. A coadded image is obtained as the average 
over multiple telescope pointings. At this stage the individual snapshots have been already 
processed with the default pipeline (versions 13.0 and upgrades), that among other basic 
operations includes the linearization and the fitting of the signal ramp, the dark subtraction 
and the fiat-fielding. Since each observation in staring mode consists of two exposures with 
different position of the source along the slit, we have subtracted the background emission 
by taking the difference between the couple of images in the nodding cycle. The spectra have 
been extracted following the standard steps for point-like sources within the software SPICE. 
The flux uncertainties have been computed considering the typical poissonian distribution 
for source and background counts. The latter are estimated from the companion coadded 
images, since the off-source subslit provides the background of the nodded on-source obser- 
vation. 

While the wavelength calibration is not a crucial issue in this work, a precise flux calibration 
is required to derive a reliable estimate of the relative AGN/SB contribution to the bolomet- 
ric luminosity (see Section [6]). In the older versions of the processing pipeline the accuracy 
of the absolute flux calibration for the Spitzer-IRS low-resolution orders was quoted to be 
about 20 per cent. Even with such a coarse evaluation our final results are not substantially 
modified, hence the Spitzer calibration as been adopted. We have anyway checked its reli- 
ability against the photometric data provided by IRAS. For 53 sources (corresponding to 
75 per cent of our sample) the 12 fim IRAS flux density is only an upper limit, which is 
always higher than the Spitzer-IKS equivalent. For 16 sources not even the IRAS flux at 
25 yum is available, while in all the remaining cases the Spitzer measure turns out to be fully 
consistent with the IRAS one (it is indeed slightly lower on average, by ~20 per cent at 
most). Even if the analysis is restricted to a narrow spectral range and only the Short-Low 
(SL) orders are involved, we performed the extraction of the Long-Low (LL) orders as well, 
so that we have the entire ~5-35 /im spectra of more than 60 sources. Only in a few cases 
a slight rescaling is necessary for a smooth connection between SLl and LL2 orders, and 
this can not introduce any systematic effect on the measured 25 fim flux. Aperture losses at 
long wavelengths can indeed occur when comparing the narrow-slit IRS spectroscopy to the 
large-aperture IRAS photometry. Whatever the right explanation may be, we conclude that 
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the Spitzer-IKS calibration at 5-8 /iin is fully reliable, and that the normalization of our 
spectra to match the IRAS fluxes, however not possible for the whole sample, is not even 
required. 



6 AGN/SB BOLOMETRIC CONTRIBUTION 

The use of spectral templates to model the intrinsic AGN and SB components significantly 
reduces the number of degrees of freedom in our analytical description. The decomposition 
method has already been introduced and discussed in detail in Paper I, nevertheless it is 
useful to examine again the main steps. By defining the 6 /im-normalized AGN and SB 
templates m^^" and m***, the observed 5-8 /im ULIRG emission can be parameterized as 
follows: 

f^bs^X) = ft [(1 - a,)ut' + ae^re-^^'^] , (2) 

where is the AGN contribution to the intrinsic (de- absorbed) flux density /g"*. Along 
with the optical depth to the active nucleus Tq = r(6 /im) and apart from the flux normal- 
ization, ttg is the only degree of freedom in our model. The best fits have been found via 
the minimization of the x^- spite of the large differences observed in our ULIRG spectra, 
this simple model always provides a good match to the data. Though the best fits have 
usually a reduced ^2 and are not formally acceptable in a statistical sense (due to the 
use of templates that ignores minor features, and the small Spitzer-IKS error bars) both 
the aromatic emission and the continuum shape are well reproduced, with residuals smaller 
than 10 per cent at all wavelengths. 

The fitting procedure allows the determination of /g"* as well, hence we can compute the 
ratio between the absorption-corrected 6 fim luminosity and the bolometric luminosity of 
each source, defined as 

--(^). (3) 

where Fjji is the total IR flux introduced in Eq.([T]). Since the equivalent quantities for pure 
AGN and pure starbursts (hereafter indicated respectively as R°-9n Q^^id R^'') are widely dif- 
ferent from each other, R is itself an indicator of the significance of AGN activity within com- 
posite sources, and can be used both to test the consistency of our decomposition method and 
to assess the relative AGN/SB contribution in terms of bolometric luminosity. By rendering 
explicit the AGN and SB components within Eq.([3]) respectively as ae/e"* and (1 — ae)/"*) 
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and decomposing as F^l"' + Ff^, the dependence of R on a^, i?"?" ^n^j ji^b ^^^^^ brought 



R = TT^. 4 



out through a simple algebra: 

R'^9^ R^^ 
a^R"^ + (1 - ae)^''^" 

This theoretical R-a^ relation has been fitted to our data, with i?'^^'" and R'^^ as fioating 
variables. In Paper I we obtained that R"'^'^ / R'^^ ~28; we have anyway decided to repeat the 
evaluation of these critical parameters, because some fits have been improved with respect to 
the previous work and three extra sources have been added to the sample. The new results 
are in good agreement with the previous ones: 

log i?"^" = -0.551°:°^ and log i?^*- = -1.91+°;°^. 

Taking R"-^"'/R'^'' ~23, we are able to provide a quantitative estimate of the AGN contribu- 
tion {ahoi = /Fm) to the bolometric luminosity of each source, as 

a^ + iR-a'-lR'^l-a^y ^' 
The values of aboi are listed in TablHl Incidentally, it should be noted that a value aboi ~ 1 
is not reached even in optically bright quasars, since the contamination from star formation 
activity is always present and non-negligible. In FigJl^ we show the relation between R and 
a^oh which can be expressed in the neat and manifest form R = atoiR"'^"' + (1 — Ciboi)R'^'' by 
inverting Eql5l 

The regularity of the pattern shaped by the location of the sources in this plot is in itself 
a qualitative validation of our method. We have anyway performed a further test of self- 
consistency, by computing for each source the following quantities: 



int ' 



\ OibolFiR J 



uq{1 - ae)/^ 



int 
6 



(^) 

Although this argument may appear as fully circular, it indeed provides a strong confirma- 
tion of our approach. In fact, as shown in FigJ3)D, the ratios -R'*^ for the SB component in 
composite sources (which depend on our AGN/SB decomposition) agree within the disper- 
sion with those for SB-dominated ULIRGs (which conversely are directly measured). Such 
agreement would be lost both in case of missed AGN detections and in case of fake AGN 
detections, since the global IR fiux is sensitive to the actual AGN contribution. 
The present method is then reliable in detecting genuine AGN components and in extrapo- 
lating the average AGN/SB contribution to the bolometric luminosities. Moreover, in spite 
of its simplicity, it proves to be effective when compared to other diagnostics at different 
wavelengths, as discussed in the next Section. It is important however to underline that the 
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AGN bolometric contribution (ex^^j) 

Figure 4. (a) Ratio R between intrinsic 6 fim and bolometric luminosity versus the AGN bolometric contribution djoi- The 
error bars of R are mainly due to the uncertainties in the total IR flux Fjn. The red solid curve traces the best fit of the R-ae 
relation from Eq|4] It is worth noting that the sources follow a very regular pattern and no evidence of outliers is found, (b) 
Same as above, with the ratios for the AGN and SB components plotted separately. The red dashed lines mark the best values 
of ij^s" and R"'' . The agreement among the SB components at different values of aj^; proves the reliability of our method (see 
the discussion in the text). 



scatter in R""^^ and R^'' is significantly larger than the uncertainty on the best values of R°'3n 
and i?'*^. The reduced of the R-a^ relation varies slowly within a region of the parameter 
space that roughly corresponds to the dispersion around the best fit in FigJl^. We then 
emphasize that such a dispersion (~0.3 dex, nearly constant with respect to ctjoi) is to be 
considered the real uncertainty in the 6 /xm to bolometric ratios for the AGN and SB com- 
ponents of individual objects (as opposed to the ensemble averages). This can slightly affect 
the estimates of aboi, increasing the statistical errors reported in TabOJas well. Nevertheless 
the results are precise enough to establish the magnitude of the AGN/SB contribution both 
to the luminosity of individual sources and to the overall energy output of local ULIRGs. 
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Table 3. Spectral parameters and AGN contribution for the 71 ULIRGs in our sample. All the optical depths refer to the AGN 
component: a colon indicates the upper limits evaluated from the residuals to the best fit. ae: AGN contribution to the intrinsic 
(i.e. de-absorbed) continuum emission at 6 fira (in per cent), ts: Optical depth to the AGN continuum at 6 fira according to the 
power-law prescription for reddening, a (,£,;: AGN contribution to the bolometric luminosity (in per cent). These entries are only 
affected by the statistical uncertainty both in the flux amplitude of the AGN/SB components and in the average ratios _R"9" 
and ij"*. The systematic effects are discussed and quantified in the text. Ti^e- Optical depth of the ~6.0 /im absorption feature 
due to water ice. Tf^ac,!-. Thac.2- Optical depths of the companion features attributed to hydrogenated amorphous carbons, 
centred at ~6.85 and ~7.25 fira, respectively. Qe*, re*, ctbol*'- Same as ae, tq and ojoj, obtained by assuming a spectral slope 
of 0.7 (instead of 1.5) for the AGN template (see Section [S]!. V/L/X: Spectral classification in the visible, L-band and hard 
X-rays. *: Starburst; ®: LINER; •: AGN detection; 0: AGN detection, tentative. 



Source 


OtQ 




T-6 








^hac,2 






re* 




V/L/X 


ARP 


220 


84il 


\ 


.47i0.01 


18+* 


< 2. 13 


90+0 01 


< 0.21: 


90itl 


2 


nn-i-n m 


20+* 

2U_3 




IRAS 


00091-0738 


97il 




.oO±U.U4 


60+^ 


n c -Ln in 
U.SoitU. iU 


U.o / itU.Uo 


0.o4±U.U4 


98itl 




. oUztU.Uo 


62it6 


^1 


IRAS 


00188-0856" 


93±1 





.37±0.04 


37 ^ 

— 5 


< 1.18 


< 0.67 


< 0.66 


93ibl 





.42 + 0.01 


29^4 


03 w — 


IRAS 


00456-2904 


< 1.0 






< 0.05 








< 1.0 






< 0.03 


,1 

— — 


IRAS 


00482-2721'' 


< 54 




< 0.04 


< 4.8 


< 0.01 


< 0.96: 


< 0. 79: 


< 38 




< 0.07 


< 1.7 


— — 


IRAS 


01003-2238 


96±1 


1 


.58±0.02 


53+? 
-5 


0.18±0.04 


0.14±0.01 


< 0.01 


98±1 


2 


.12+0.02 


56+^ 


©1 - - 


IRAS 


01166-0844 


> 99 


2 


„.+0.23 
■■^=-0.08 


89+'^ 


41 + 0-36 
^■^-^-0.33 


70+"-3S 
'1-0.35 


40+''-l^ 
"-*"-0.18 


> 99 


2 


qo-i-n HA 


90+3 


y,l 


IRAS 


01298-0744 


> 98 




'7Ci-i-r\ no 


77+5 


U.yiztU.Uo 


A .1 7_l_n riQ 
U.4 / itU.Uo 


n on^n no 


> 99 




OH-Ln no 


78+5 


^1 


IRAS 


01569-2939 


85±1 


1 


1 O + 002 
■"Io.04 


13-3 


n -50 + 0.09 


0.31±0.09 


0.31±0.05 


91±1 


1 


.66 + 0.05 


22^4 


,1 

ik — — 


IRAS 


02021-2103 


78±2 


1 


.91±0.10 


13+'l 
-3 


< 0.01 


< 0.40: 


< 0.18: 


86±1 


2 


.47+0.11 


15+* 
-3 




IRAS 


02411-1-0353= 


< 17 






< 0.9 








30±2 





.77+0.07 


1 + 0.4 
1-2-0.3 


,1 

— — 


IRAS 


03250-1-1606 


< 3.4 




< 0.11 


< 0.2 








< 3.9 




< 0.14 


< 0.2 


/ry.1 ^4 
— 


IRAS 


04103-2838 


59±1 





.14±0.02 


5-9-1.1 


< 0.01 


< 0.01 


< 0.01 


70itl 





.63+0.02 


''-3-1.1 


rTil .9 


IRAS 


05189-2524 


92+1 




< 0.01 


Lfl 

32tio 


< 0.15 


0.04±0.01 


< 0.01 


94±1 





.43+0.01 


30^= 


14 10 


IRAS 


u 1 oyo-f-Douo 


> 99 




< 0.06 


86±4 


< 0.01 


< 0.01 


< 0.02 


> 97 




< 0.01 


94+* 


,2 ,4 Qll 


IRAS 


08559-1-1053 


68±1 




< 0.01 




< 0.01 


< 0.16 


< 0.12: 






< 0.01 


5 0+1-2 
^-"-2.4 


.1 .4 - 




uoD ( ^-f-oyio 


> 99 





.44±0.01 


87±3 


0.16±0.02 


0.24±0.01 


O.lOiO.Ol 


> 99 





.91+0.01 


81±4 


0l - 


IRAS 


09039-1-0503 


61±1 





.71±0.04 


6 4+1'* 
" -1.3 


< 1.33 


< 0.01 


< 0.38: 


72±1 


1 


.20+0.04 


6 g+i-** 


0l 0" - 


IRAS 


09116-1-0334 


< 1.8 




- 


< 0.08 


- 


- 


- 


< 0.9 




- 


< 0.03 


0l 04 - 


IRAS 


09539-1-0857 


91±1 


1 


.85±0.03 




< 1.27 


0.31±0.03 


0.28±0.02 


94±1 


2 


.38+0.03 


33^5 


0l - 


IRAS 


10190-1-1322 


< 0.3 




- 


< 0.02 


- 


- 


- 


< 0.3 




- 


< 0.01 


*1 i,"- - 


IRAS 


10378-1-1109 


80±1 





.21±0.02 


15^3 


< 2.56 


0.43±0.03 


0.48±0.03 


87±1 





.70+0.02 


16+3 


0l 0'' - 


IRAS 


10485-1447 


60il 





. 16i0.03 


•=■2-1.2 


< 2.00 


A CA_I_A A7 

O.oOihO.O I 


< 0.34: 


70itl 





.61±0.04 


''-2-1.2 


— 


IRAS 


10494-1-4424 


< 0.4 






< 0.02 








< 0.4 






< 0.01 


0i 0'' - 


IRAS 


11095-0238 


97±1 


1 


■ 


63±5 


• i i 






98itl 


2 


■"^-0.04 


KO+11 
''2-8 


^-7\l j,4 


IRAS 


11130-2659 


84±1 


1 


.15±0.02 


1 Q + 4 


0.37±0.05 


0.70it0.05 


0.40±0.07 


90itl 


1 


.67±0.02 


21+4 


0l 


IRAS 


11387-1-4116 


< 0.8 






< 0.04 








< 0.7 






< 0.02 


^1 ,4 


IRAS 


11506-1-1331 






< 0.01 




< 0.29 


n Qo + 0.02 
0-32-0.20 


< 0.01 


57itl 





.18±0.02 


3-6-0.7 


-*- — 


IRAS 


12072-0444 


95±1 


1 


.07±0.01 


.,0+6 
43^" 

— 5 


0.28±0.02 


0.26±0.01 


0.19±0.01 


97±1 


1 


.58±0.01 


45^ 


.1 .4 
• • — 


IRAS 


12112-1-0305 


< 13 




< 0.07 


< 0.7 








< 9.7 




< 0.04 


< 0.4 


^1 5 12 


IRAS 


12127-1412= 


> 99 




< 0.04 


S9+1 


< 0.44 


< 0.40 


< 0.38 


> 99 





.37±0.01 


8'4_4 


^-7\1 .4 
09 • — 


IRAS 


12359-0725 


54 

-43 




< 0.01 


A n + 1 6 

^^-4.4 


< 0.37 


< 0.48 


< 0.38 


ro+1 

53-8 




< 0.01 


1 +0 s 
3-1-1.2 


0l 0'! — 


IRAS 


13335-2612 


< 0.6 






< 0.03 








< 0.5 






< 0.02 


^1 

09 


IRAS 


13454-2956 


61 + 1 
" -27 




< 0.01 




< 0.01 


< 0.16: 


< 0.09: 


59±1 





.07±0.04 


4 0+1-1 
"-0.S 


.1 - - 


IRAS 


13509-1-0442 


< 0.6 






< 0.03 








< 0.6 






< 0.02 


*1 


IRAS 


13539-1-2920 


< 0.4 






< 0.02 








< 0.4 






< 0.02 


*1 k^- 


IRAS 


14060-1-2919 


< 0.4 






< 0.02 








< 0.4 






< 0.01 


*1 k*- 


IRAS 


14197-I-0813-' 


75±2 


2 


.10±0.14 


12;* 


< 1.61 


< 0.01 


< 0.01 


84±2 


2 


.63±0.15 


13;* 




IRAS 


14252-1550 


< 20 




< 0.04 


< 1.1 








< 24 




< 0.31 


< 0.9 


0l *4_ 


IRAS 


14348-1447 


51+1 




< 0.09 


43+1.7 


< 1.94 


"-^'-0.36 


62+''-l** 
"-"^-0.26 


59±3 





.25±0.10 


3 9+1-5 
^-'^-l.o 


0l *5 ^12 


IRAS 


15130-1958 


91 + 1 

^ -3 




< 0.01 


on + 6 
3^-10 


< 0.08 


< 0.01 


< 0.01 


94±1 





.42±0.01 


30^^ 


.1 .4 - 


IRAS 


15206-1-3342 


52±1 





.30±0.02 


4 5+1-2 


< 0.01 


< 0.01 


< 0.01 


65±1 





.83±0.02 


5 1+1-3 


*1 


IRAS 


15225-1-2350 


87±1 





.74±0.02 




0.92±0.05 


0.35±0.02 


0.22±0.01 


92±1 


1 


.26±0.02 




*1 04- 


IRAS 


15250-1-3609 


96±1 


1 


.06±0.01 


53±6 


1.38±0.03 


0.60±0.03 


0.41±0.02 


98±1 


1 


.58±0.01 


55+8 


02 _ *12 


IRAS 


15462-0450 


90+1 




< 0.01 




< 0.01 


< 0.07: 


09+''-°l 
"-"'^-o.os 


89±1 





.22±0.01 


19+* 


.1.4- 


IRAS 


16090-0139 


89±1 





.69±0.01 


26±4 


0.90±0.03 


0.57±0.02 


0.36±0.01 


93±1 


1 


.21±0.01 


28+= 


0l 04 - 


IRAS 


16156-1-0146 


94±1 





.61±0.02 


43+1 


0.48±0.04 


0.24±0.01 


O.lliO.Ol 


96±1 


1 


.09±0.01 


43;^ 


.1 - - 


IRAS 


16468-1-5200 


85±1 





.77±0.02 


20tt 


< 1.50 


0.94±0.08 


0.60±0.04 


91±1 


1 


.30±0.02 


22tl 


0l*4_ 


IRAS 


16474-1-3430 


< 4.9 






< 0.3 








< 5.5 




< 0.01 


< 0.2 


*1 04- 


IRAS 


16487-I-5447S 


21^1 




< 0.04 


1.2±0.4 








31±2 





.47±0.08 


1 2+°-'' 
1-2-0.3 


0l 04 - 


IRAS 


17028-1-5817 


< 1.2 






< 0.06 








< 1.1 






< 0.04 


0l 04 - 


IRAS 


17044-1-6720 


91±1 





.32±0.01 


29+4 


0.09±0.02 


< 0.01 


< 0.01 


94±1 





.82±0.01 


3lt= 


0l .4 - 


IRAS 


17179-1-5444 


84±1 





.31±0.02 


18;* 


< 0.01 


< 0.01 


< 0.01 


89±1 





.81±0.02 


19;* 


.1.4- 


IRAS 


17208-0014 


< 7.9 






< 0.4 








< 7.5 






< 0.3 


03 *5 ^12 


IRAS 


19254-7245'' 


89±1 





.21±0.08 


26+« 

— 5 


21 + ''-l" 


0.21±0.04 


0.10±0.02 


88±1 





.09±0.01 


17+3 


.3 .5 ,12 



20 E. Nardini et al. 

Table 3 - continued Spectral parameters and AGN contribution for the 71 ULIRGs in our sample. 
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References. ^ Veilleux et al. (1999a), ^ Veilleux et al. (1995), ^ Due et al. (1997), * Imanishi, Dudley & Maloney (2006), 
5 Risaliti et al. (2006b), ^ Risaliti et al. (2006a), ^ Imanishi, Dudley & Maloney (2001), * Iwasawa et al. (2005), ^ Teng et 
al. (2008), 1° Severgnini et al. (2001), " Imanishi & Terashima (2004), Franceschini et al. (2003), Balestra et al. (2005), 
14 Vignati et al. (1999), Imanishi et al. (2003). 

Notes. " The value of tq is obtained allowing Tagn = 0.78. This is indeed considered an AGN detection. A comparable 
minimum in the parameter space suggests a possible AGN component with re = 0.24 and ctjo; = 1.1 (per cent). The value 
of re is obtained allowing Fagn = 0.58. The value of tq is obtained allowing Fagn = 1.31. ^ A comparable minimum in the 
parameter space suggests that this source could also be SB-dominated. ^ Due to the dispersion around the templates, this 
value of ojoi is too small to include this source among the safe AGN detections. The value of rg is obtained allowing 
^agn = 0.52. * The value of rg is obtained allowing Tagn = 1.02. 



In the light of the elements collected so far, we are able to address again the question 
about the existence of extremely obscured AGN components that can not be detected even at 
5-8 /im, but only contribute to the IR emission at longer wavelengths (e.g. at ~25-60 fim). 
Such components would pull downward our estimate of logi?'*'' = —1.91. An empirical 
measure of this quantity can be performed at lower IR luminosities on the B06 starburst 
sample, yielding an average of logW'' = —1.66. This latter value anyway requires some 
caution for two opposite reasons: some objects in B06 (e.g. NGC 1365, NGC 4945) show 
significant nuclear activity; on the other hand, we are dealing with nearby sources and 
aperture losses in the SL spectral orders can occujf]. Disregarding such complications (these 
systematic effects are believed to compensate for each other), it turns out that the two 
distributions of i?*^ fairly overlap, due to their large dispersion. Even if they appear to be 
somewhat displaced with respect to each other, the variation is too small to support the 
hypothesis of a completely buried AGN population escaping this 5-8 /xm spectral probe (or 



^ We have adopted the same strategy of B06, scaling up all the spectral orders to match the LLl flux density, and then 
normalizing to IRAS at 25 /im. By comparing the final spectra to the 12 fj,m IRAS fluxes we evince that a residual loss may 
be present, and provide an upper limit for the average of logi?'' « —1.5. 
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any other Spitzer-IRS probe; see also Veilleux et al. 2009). This difference can be possibly 
ascribed to the overall physical properties of the starburst process itself; for example, in the 
ULIRG luminosity range the starburst is triggered by a major merger (e.g. Dasyra et al. 
2006), while at lower luminosities it mostly takes place in isolated systems. 

7 COMPARISON WITH OTHER DIAGNOSTICS 

Our 5-8 fim AGN/SB decomposition has provided 50 convincing AGN detections out of 71 
sources in our sample, confirming all the already known AGN components and uncovering 
some more. Conversely, the AGN detection rate among ULIRGs is ~30 per cent at optical 
wavelengths (Veilleux et al. 1999a), and < 50 per cent according to mid-IR high-ionization 
lines (Farrah et al. 2007). X-ray diagnostics is actually extremely powerful, but can be 
applied on such a large scale only collecting a huge integration time: after Teng et al. 2005, 
we estimate that an average exposure of ~30 ks per source is required in order to reach up to 
our detection rate with the present X-ray facilities. In Section [2] we have briefly summarized 
the multiwavelength picture of ULIRGs which is obtained from the most effective diagnostic 
methods used so far. We now discuss in more detail how our results compare with those of 
previous works, allowing us also to overcome some of their limitations. 

7.1 Optical 

The nuclear regions of ULIRGs are affected by typical obscurations of several tens of mag- 
nitude at visible wavelengths, hence it is difficult to interpret the results of the diagnostics 
based on emission lines. None the less some ULIRGs exhibit in their spectra the broad 
optical recombination lines characteristic of Seyfert 1 galaxies. Even when a direct access 
to the broad line region is obstructed by larger amounts of dust, possibly aggregated in a 
toroidal shape as expected for type 2 Seyfert-like objects, the presence of a working AGN 
can be firmly established in the near-IR via the detection of broad permitted features or 
of the high-ionization 1.962 /xm [Si vi] line (Veilleux, Sanders & Kim 1999). In principle, 
a hidden broad line region can also be looked for in polarized light, a technique that has 
been developed with success for nearby Seyfert 2 galaxies (after Antonucci & Miller 1985) 
but is almost unavailable for ULIRGs (e.g. Pernechele et al. 2003) because of their larger 
distance, and possibly frustrated by the greater complexity in the geometrical structure of 
the dust (large covering factors and dumpiness are not easily reconcilable with the reflection 
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Figure 5. Optical depth to the AGN component evaluated at 6 fim versus the AGN intrinsic contribution to the 5-8 fim 
emission. Each source is plotted with a different symbol according to its optical classification, as indicated in the box. The 
vertical black dotted lines mark an AGN contribution to the bolometric luminosity of 1, 10 and 25 per cent, respectively. The 
crosswise orange curve represents our empirical detection limit: due to the dispersion around our AGN/SB templates, an AGN 
component on the left of this line could be missed by our method (because its observed contribution is too faint). For the same 
reason, AGN detections close to this region of the plot are to be considered as tentative. The bolometric contribution of such 
components would be anyway negligible (a few per cent at most). 



scenario, nor with dichroic transmission). In most cases only narrow features are observed, 
and it is necessary to resort to line ratios in order to classify a ULIRG. The separation 
of active nuclei and star-forming galaxies is based on empirical boundaries, discriminating 
the effects of photoionization by hot young stars from those due to non-thermal sources 
with a power-law continuum. This is in itself a limitation to the accuracy of this diagnos- 
tic method, which is further degraded by the extinction effects. According to Veilleux & 
Osterbrock (1987), three equivalent diagnostic diagrams can be set up using the line ratios 
[O l] A6300/Ha, [N ll] A6583/Ha and [S ll] (A6716+A6731)/Ha along with [O ill] A5007/H;g 
(the latter distinguishes between a type 2 Seyfert galaxy and a LINER). It turns out that a 
significant number of ULIRGs do not belong to the same spectral type in all these diagrams, 
with some degree of overlap among LINERs and H II regions (Veilleux et al. 1999a). In any 
case, no evidence of AGN activity is clearly found within these two classes. 
In Figl5] the AGN intrinsic contribution to the 5-8 fim emission and its obscuration, as 
obtained with our technique, are plotted together with the optical properties. With some 
scatter, the location of each source correlates with its optical classification. In particular, 
there is a tight agreement for Seyfert-like ULIRGs: type 1 objects are found in the bottom 
right-hand corner, with a large AGN content and vanishing obscuration, while the average 
type 2 hosts a slightly dimmer and obscured AGN. The interpretation of the plot is less 
straightforward when one considers also the sources that are classified as H ll regions or 
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LINERs. It is clear that optical spectroscopy, besides giving incomplete information about 
the relative AGN/SB contribution, can also be misleading: this is indicated by the presence 
of only one type 2 Seyfert-like source among the nine entries with aboi > 0.25 and rg > 1. 
Also four H II galaxies fall into this group, and they are all true AGN detections at a sim- 
ple visual inspection, sharing a reddened and absorbed continuum and strongly suppressed 
aromatic features. These targets are hence perfect candidates for follow-up observations at 
X-ray wavelengths. As opposed to the unexpected scatter of H ll galaxies, the wide disper- 
sion of LINERs in FigjS] confirms again that such systems are rather heterogeneous with 
respect to the nature of their energy source. Out of 33 LINERs comprised in our sample, 23 
contain an AGN component, with a wide range of bolometric contributions. This fraction 
(~70 per cent) is very close to the detection rate achieved by collecting all multiwavelength 
diagnostics, but is significantly higher than those obtained from radio, optical and X-ray 
observations separately. Our decomposition method is then capable to resolve most of the 
ambiguity concerning the ionization mechanism in LINERs (see also Dudik, Satyapal & 
Marcu 2009; Gonzalez-Martin et al. 2009). 

7.2 Mid-IR 

A wealth of studies at mid-IR wavelengths have been published since the launch of the In- 
frared Space Observatory (ISO; Kessler et al. 1996), most of which rely on different mixtures 
of three basic diagnostics, i.e. hot dust continuum, aromatic emission and silicate absorption. 
The usual strategy makes use of single indicators such as the equivalent width of the 6.2 
or 7.7 /im PAH feature, flux ratios involving fine-structure lines, the 9.7 fim silicate optical 
depth, the 15 to 6 /im continuum colour. Pairs of these diagnostics are used to construct 
bidimensional diagrams in which the location of ULIRGs is compared with those of stan- 
dard AGN and SB control samples (e.g. Genzel et al. 1998). Alternatively, one can attempt 
at reproducing the whole observed emission within a certain wavelength range by means 
of a limited number of components. Both spline continua plus variable PAH emission and 
fixed spectral templates have been employed. Our method follows the latter approach, and 
is qualitatively similar to those presented in Laurent et al. (2000) and Tran et al. (2001), 
with the great improvement in data quality made available by Spitzer. It is difficult to find 
a quantitative correspondence among the results obtained on the individual objects, since 
these earlier works were intended primarily to derive the statistical properties of the ULIRG 
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Figure 6. Equivalent widths of the 6.2 /im (blue diamonds) and the 7.7 fim (red rings) PAH features (from Veilleux et al. 2009) 
plotted against our estimate of the AGN observed contribution to the 5-8 fim continuum, evaluated at the same wavelength. 
The typical error bars are comparable to the symbol size. 

population as a whole. A possible example of a more specific comparison is illustrated in 
FigEl where the AGN observed contribution to the 5-8 fim emission of 48 sources is com- 
pared to the equivalent widths of the 6.2 and 7.7 fim aromatic features (from Veilleux et 
al. 2009). The expected anti-correlation between these quantities can be appreciated, and 
the scatter is highly reduced with respect to the measurements performed before the Spitzer 
era (e.g. Rigopoulou et al. 1999). At the same time it is clear that reddening extensively 
shapes the spectra of ULIRGs, and its effects have a substantial impact when estimating 
the AGN/SB bolometric contribution. 

Another way to probe the degree of absorption and the geometrical structure of the dust, 
and test the characterization we offer of these issues, is provided by the silicate feature at 
9.7 fim: its opacity, together with the intensity of aromatic emission, can be plotted in a dia- 
gram useful not only for working as a classification scheme, but also for tracing the possible 
evolutionary path of a source. The silicate feature has not been studied in the present work. 
We anyway point out the good agreement between our findings about the most notable 
sources in our sample and their position in such a diagnostic diagram (Spoon et al. 2007). 
In particular, the silicate strength appears to have a positive correlation with the reddening 
we have estimated for the AGN component. 

7.3 High-ionization lines 

As mentioned before in relation with the 1.962 fim [Si vi], the presence of fine-structure 
lines from highly-ionized atoms has a powerful diagnostic value in distinguishing between 
black hole accretion and star formation in ULIRGs. Also in the mid-IR, the mere detection 
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of features such as the [Ne v] A14.32 hue is the clear mark of an active nucleus, because of 
the required ionization energy of 97.1 eV. Again, by combining suitable line ratios among 
themselves or with different AGN/SB indicators, useful diagnostic diagrams can be set up, 
that make it possible to point out the main average trends but can be misleading if applied 
to single objects. 

Recently Farrah et al. (2007) have analyzed the high-resolution Spitzer-IRS spectra of 53 
local ULIRGs, 27 of which are also included in our sample. This allows a meaningful compar- 
ison between our AGN/SB decomposition method and diagnostics based on fine-structure 
lines. Focusing on the common sources, three of them have been used for the construction of 
our SB template, while signatures of AGN activity have been found in each of the remaining 

24 according both to Paper I and the present work. We therefore look for a confirmation of 
this activity through high-ionization lines. The simplest approach is of course the detection 
of individual features. [Ne v] A14.32 is found in ten objects, in nine of which we detect an 
AGN component as well. The only exception is IRAS 20414—1651, that we have assumed to 
be a prototypal starburst because of its prominent PAH features and cold mid-IR colours. 
IRAS 20414—1651 is actually the only case without the simultaneous detection of [O iv] 
A25.89, another usual clue of hard radiation fields. On the contrary, IRAS 10378+1109 ten- 
tatively shows the oxygen line but no neon emissioio- It is not obvious to relate the absence 
of [Ne v] A14.32 to a higher degree of obscuration. Both MRK 231 and IRAS 15462-0450, 
for instance, are not detected in [Ne v] A14.32 despite their nature of type 1 Seyfert galaxies 
at optical wavelengths. The same holds for IRAS 07598+6508. None the less, it is worth not- 
ing that the median value of the optical depth to the AGN component is ~0.3 for the nine 
ULIRGs with a [Ne v] A 14. 32 detection and ~0.5 for the remaining 15 composite sources, 
and indeed the argument of a larger obscuration seems to be reasonable for many of them. 
A high covering factor of the absorbing dust can even prevent the formation of a narrow 
line region and the production of fine-structure lines. A complete revision of both line fluxes 
and upper limits can be found in the recent work of Veilleux et al. (2009), that also includes 

25 additional sources within our sample. The considerations made so far are widely corrob- 
orated. There are six more objects with a [Ne v] A14.32 line: only one (IRAS 13335—2612) 
is classified as SB-dominated after our 5-8 /xm spectral analysis. Also the loose correlation 



^ This detection is not actually confirmed by Veilleux et al. (2009), as well as the detections in 4C +12.50; the [Nc v] A24.32 
fine is instead found in IRAS 15206+3342. 
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between successful detections and little opacities is possibly emphasized. 
The line ratios allow more complete diagnostics in combination with some other AGN/SB 
indicator over the electromagnetic spectrum. All the diagrams proposed by Farrah et al. 
(2007) point to star formation as the dominant power supply in ULIRGs: the active nucleus 
is as powerful as the starburst in only ~20 per cent of the sources, and of lesser significance 
in more than half. These fractions are in perfect agreement with our estimates. The effective- 
ness of mid-IR fine-structure lines in ULIRG diagnostics is anyway restricted to the finding 
of general trends, since no firm quantitative constraint can be put on the AGN bolometric 
contribution (except for some remarkable examples), and even qualitative considerations 
may lead to mistaken conclusions on individual objects. 

7.4 Hard X-rays 

Since ULIRGs are faint X-ray emitters, with a 2-10 keV continuum fiux ranging from ~10~^ 
to a few xlO~^ of their overall IR output, only a limited number of them has been observed 
successfully with the Chandra and XMM-Newton sateUites. We can therefore provide only a 
brief compilation of the main results that have been obtained so far. Clear X-ray signatures 
of AGN activity are a 2-10 keV luminosity exceeding 10^^ erg and a large column den- 
sity {Nh > 10^^ cm~^) absorbing an unresolved component. The detection of the iron Ka 
fluorescent line at 6.4 keV with an equivalent width of ~1 keV also indicates the presence of 
an obscured active nucleus, whose primary radiation is reflected by a cold mirror of gaseous 
material (for a detailed review about this subject see Reynolds & Nowak 2003). All these 
indicators have been found in IRAS 19254-7245 (Braito et al. 2003, 2009). MRK 231 lacks 
a strong iron line, but the direct AGN emission emerges above 20 keV (Braito et al. 2004); 
moreover this source exhibits significant nuclear variability on timescales of a few hours 
(Gallagher et al. 2002). Both nuclei of NGC 6240 have hard X-ray emission that possibly 
dominates the global energetics, with heavy absorption and a prominent neutral iron line 
(Vignati et al. 1999; Komossa et al. 2003). IRAS 05189-2524 is a typical Compton-thin 
Seyfert 2 galaxy, in which the estimated luminosity of the AGN component can not ac- 
count for the huge IR emission (Severgnini et al. 2001); also MRK 273 (Balestra et al. 2005) 
and UGC 5101 (Imanishi et al. 2003) are likely SB-dominated, although an AGN is detected 
above ~3-4 keV. The hard X-ray luminosity and spectral properties suggest an AGN contri- 
bution to IRAS 20551-4250 and, tentatively, to IRAS 23128-5919 and IRAS 20100-4156, 
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while no convincing evidence of AGN is found in IRAS 12112+0305, IRAS 14348-1447, 
IRAS 15250+3609, IRAS 17208-0014 and IRAS 22491-1808 (Franceschini et al. 2003; 
Ptak et al. 2003). 

Even when measured, the luminosity of the above-mentioned AGN components generally 
lies in the range of Seyfert galaxies. In order to represent the major energy source of the 
ULIRG bolometric luminosity they should be located behind a column density of obscuring 
material approaching 10^^ cm~^. It is still matter of debate if this is the case for ARP 220, 
the nearest ULIRG, inside which an elusive active nucleus is suspected to be at work. A 
strong iron hue has indeed been detected in ARP 220 by Iwasawa et al. (2005), challeng- 
ing the origin of the 2-10 keV emission from X-ray binaries and hence the pure starburst 
interpretation (for the X-ray emission processes in starburst galaxies see Persic & Rephaeli 
2002). 

A handful of additional sources has been recently studied by Sani, Risaliti & Salvati (in 
preparation) using also the BeppoSAX observations up to 100 keV. The most striking re- 
sult concerns IRAS 12072—0444, which is detected only at energies beyond 20 keV and 
shows up as a Compton-thick source harbouring a type 2 quasar with a column density 
Nh ~4x102^ cm-2. IRAS 08572+3915 is not detected but is suggested to be a type 2 
quasar as well, whose reflection efficiency is very low by virtue of a cocoon-like obscuring 
geometry. Finally, on the basis of the ratio between the 2-10 keV and the IR emission, 
IRAS 10190+1322 and IRAS 10494+4424 confirm their SB-dominated nature, while the 
presence of an AGN component can not be ruled out in IRAS 00188-0856, IRAS 01003-2238 
and IRAS 16090—0139: according to their Spitzer-IRS spectra and the findings of the present 
work these latter objects are believed to host a buried AGN, and are expected to be ab- 
sorbed even at 10 keV. Regrettably one of them (IRAS 16090—0139) was observed but not 
detected by BeppoSAX, and the upper limit to its flux above 20 keV is not constraining; the 
other two have not been observed at all. 

In conclusion, our results are in qualitative agreement with all the pieces of evidence from 
X-ray studies. The X-ray diagnostics of ULIRGs is very powerful, yet our 5-8 //m analysis 
provides fully consistent indications and, at present, gives access to a much larger number 
of sources. 
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8 DISCUSSION 

Having illustrated the global results of our study and their consistence with the multiwave- 
length picture of ULIRGs, in this Section we discuss in quantitative terms the reliability of 
alternative AGN templates. We also address some uncertain or anomalous cases, and their 
consequences with respect to the dust properties within the ULIRG class. We finally inves- 
tigate whether our analysis is able to recover the claimed correlation between AGN activity 
and bolometric luminosity. 

8.1 Fitting of peculiar sources 

In our model the great diversity of the observed ULIRG spectra can be entirely ascribed 
to the relative amount of the AGN contribution and its obscuration, intended as both ab- 
sorption and reddening. The main features are adequately reproduced and a good fit is 
obtained for most of the spectra in our sample, with only a few exceptions whose analy- 
sis required additional assumptions. The main issue is the possible fiattening of the AGN 
component (as found by N07), which entails a relaxation of our general template. The rep- 
resentative case is that of IRAS 07598+6508, which is unambiguously an AGN-dominated 
source: its continuum trend can be reproduced by simply allowing the AGN spectral index to 
be ~0.6. IRAS 19254—7245 (the Superantennae) is instead a composite system with regular 
absorption features. In order to obtain a good fit we need again to provide a fiatter AGN 
component: including its intrinsic slope among the degrees of freedom we obtain a value of 
~0.5, with a 6 yum optical depth of ~0.2. 

The remaining three particular cases (that in the following will also be labelled as anomalous) 
are more interesting, in that the supplementary presence of a broad and deep absorption over 
most of the 5-8 fim range causes the failure of our continuum-based AGN/SB decomposition. 
All these objects are likely to host a completely buried AGN; their emission is shown in FigJTl 
where they are compared to the benchmark spectrum of IRAS F00183— 7111 (Spoon et al. 
2004; Nandra & Iwasawa 2007). Qualitatively, IRAS 12127-1412 is the most striking object 
in our sample, since the broad absorption completely replaces the 6.2 fim PAH emission (no 
trace of which is clearly visible) and stretches longward of 7 fim. In IRAS 00188—0856 and 
UGC 5101 a strong, hard spectrum (closely resembling the hot dust continuum of AGN) 
stops abruptly: this occurs in some other sources as well, but the putative AGN continuum 
is rather steep and only a proper modelling of the absorption profiles is entailed. On the 
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Figure 7. The puzzling spectral trend of the anomalous sources in our sample, compared to the 5—8 fim emission of IRAS 
F00183— 7111 {thin red line, scaled to match the putative hot dust continuum). The apparent flatness can indeed be connected 
to extreme absorption of the AGN component within a different extinction scenario. 

contrary, the three cases at issue suggest once again a flattening of the AGN component 
with respect to rag?i=1.5, and the aromatic features, if present, seem to emerge from the 
bottom of a deep trough. A customized fitting procedure is therefore needed. The observed 
continuum is assumed to be AGN-dominated at ~5.5 fim, and since the available window 
is very narrow we have employed the 15 fim flux as a pivotal point of the fit. The AGN 
intrinsic slopes obtained in such way range from ~0.8 to ~1.3, with tq between ~0.4 and 
^0. We have then frozen the continuum and included the SB component and the absorption 
profiles in the next fitting step. The silicate feature is also required longward of ~7.5 /xm, 
as in other few objects. The physical interpretation of these anomalous spectra, concerning 
both flatness and extreme absorption, will be discussed in detail in the following subsections. 

8.2 The AGN intrinsic spectral slope 

A power-law behaviour of the 5-8 /im hot dust continuum in AGN is confirmed by the ob- 
servations, which however provide different values of the spectral index. We recall that the 
slope of the AGN template can not be treated as a free parameter, because of its degeneracy 
with the optical depth. Our assumption of a steep intrinsic gradient has been explained in 
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Section [31 Conversely, we have just presented a few peculiar cases whose spectra are better 
reproduced by a lower value of T agn- It is then worth testing how our estimates of the AGN 
bolometric contribution depend on this choice. 

Since the quahty of our fits is good and the SB template is not modified, we expect the 
observed AGN component to be nearly invariant. Hence the flatter the intrinsic continuum, 
the larger the optical depth, and consequently the AGN flux amplitude. As a first approxi- 
mation, the relative weight of the AGN component is then increased after the flattening of 
the intrinsic spectral shape. We have performed again the fitting of all our sources, modify- 
ing the AGN template and assuming the new value of Vagn=^-^ for its slope, as suggested 
by N07. The goodness of our fits remains virtually unchanged, and no convincing additional 
AGN detection is found. As forecast, however, both the optical depth and the AGN ampli- 
tude increase remarkably, because of the degeneracy between T agn and tq connected to the 
assumed recipe for the extinction law. The systematically larger AGN flux also affects our 
estimate of i?"^^", from the fitting of the R-a^ relation of Eq.(jll). The modified values of the 
6 /xm to bolometric ratios are the following: 

log^asn _ _Q_3g+0.06 \ogR'^ = -IMt^^^^. 

With respect to the estimates of Section[6l R^^ is confirmed (the SB template is not changed), 
while R"'^'^ turns out to be ~0.2 dex higher. The revised values of a^oi are again computed 
from Eq.([5]), where now i?"^"//?*^ ~35. A striking result is obtained as a consequence of 
the correlation between ag and i?*^^", i.e. a moderate variation in the slope of the AGN 
template has no significant effect on Uboi- This invariance can be better appreciated in FigJHl 
while Tabl3] lists the two alternate values of both the AGN bolometric contribution and the 
optical depth. This is of course a notable finding. Nevertheless we emphasize once again that 
our original template appears to be more suitable to describe the AGN component inside a 
ULIRG, in spite of the peculiar cases. The alternative scenario, in fact, implies that virtually 
all the detected AGN components are heavily reddened; as a consequence, the distribution 
of continuum opacities is discontinuous and strongly biased towards the highest values. This 
anomaly can not be explained even allowing a far larger dispersion in the 5-8 yum AGN 
spectral shapes than that observed by N07. 
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Figure 8. Different estimates of the AGN bolometric contribution (ojo; and Oj,,,!*) as a function of its intrinsic spectral slope. 
The vertical axis shows the values obtained for ragn=l-5, the horizontal one those obtained for ragn=0.7. The red solid line 
marks the equality between the two estimates, while the blue dashed lines indicate a deviation of 50 per cent: each entry falls 
inside this region. For clarity, 17 sources that are confirmed to be SB-dominated in both cases have been omitted, and fall 
outside the region plotted above. We recall that, due to the dispersion around our templates, safe AGN detections require 
atol >0.01. 

8.3 Absorption features 

We have already mentioned liow tlie SB component is affected by internal extinction only; 
the large spread in the optical depth of the 9.7 and 18 fim silicate absorption features within 
the SB class suggests a wide range of dust properties, a full investigation of which would 
require a larger spectral coverage and a more detailed model than those considered in this 
work. Concerning the AGN component, it is not possible to put constraints on the nature 
or geometry of the cold screen blocking the view to the nuclear region. However, we are able 
to discuss the properties of the main absorption features and the accuracy of a power-law 
extinction curve. 

In general, the AGN and SB templates we have adopted are sufficient to account for all the 
major features observed in the 5-8 fim spectra of ULIRGs. Minor features, when involved, 
have been reproduced through gaussian profiles, but do not alter the AGN/SB decomposi- 
tion. There are indeed plenty of possible absorption features falling in the 5-8 fim range, as 
can be also inferred from the observations of a highly-obscured environment like the Galac- 
tic Centre. Notably, none of these has been found in our spectra without the simultaneous 
detection of a significant AGN contribution to the continuum emission. A characteristic 
couple of features is frequently distinguished in composite sources, centred at ~6.85 and 
~7.25 and ascribed to the C-H bending modes of aliphatic hydrocarbons (also known 
as hydrogenated amorphous carbons, HAG). They usually display a regular gaussian profile, 
and their relative depth can be used to constrain the chemical properties of the dust grains 
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(see also Dartois & Munoz-Caro 2007). The most sensitive point for our study is anyway 
the possible presence of another absorption feature, centred at ~6.0 fim, that is commonly 
attributed to a mixture of ices and is much more complex. Although this profile can be 
strongly asymmetrical, with an elongated red wing, it is usually so shallow and/or regular 
that a rigorous selection of its shape is not required. In the mentioned case of Sgr A* the 
6 /xm feature would be optimally fitted by means of pure H2O ice at 30 K (Chiar et al. 
2000), and this consequently represents a reasonable modelling for our purposes. We there- 
fore decided to adopt the pure water ice profile at 30 K (provided by the Leiden Observatory 
database) as the backbone of the 6 fim absorption in our fits. 

As discussed in the following, this assumption is challenged in only a few objects hosting 
a deeply enshrouded AGN (including the three anomalous spectra, where it actually fails). 
The plain identification in terms of water ice, in fact, would lead to a large overestimate 
of the optical depth in the companion 3 fim absorption band. A simultaneous fit of both 
profiles in Sgr A* drives only a partial explanation of the 6 /im feature as the product of a 
mixture of ices, i.e. HaOiNHgiCOa (100:30:6) at 15 K and pure HCOOH at 10 K (Chiar et 
al. 2000). Most of the residuals are accounted for by a 6.2 fim component whose origin is 
supposed to be aromatic. This remarkable example proves how the exact nature of the whole 
6 fim feature is still poorly known. Yet the formation by surface chemistry of many different 
species in the icy mantles of the dust grains suggests the possibility of large variations in its 
profile, as a function of the chemical and physical properties of each source. This is indeed 
what we find out, since the central wavelength (actually the wavelength of minimum flux) 
appears to shift within the ~5.9-6.2 fim range, and also the blue and red wings tend to be 
slightly different for the objects whose optical deptlJll (hereafter r^ce) is larger than that of 
Sgr A*. A sizable value of Tjce can be considered itself a qualitative indicator of a significant 
AGN contribution, but in extreme cases it can even prevent a straightforward AGN/SB 
decomposition. 

In the light of the above considerations, a likely interpretation of the broad absorption band 
characterizing some of our 5-8 fim spectra is the blending of the 6 fim ice profile with the 
pair of HAG features around 7 fim. Evidence of a dramatic broadening of the 6.85 fim line 
is actually found in the observations of several embedded protostars (Keane et al. 2001). 



^ All the optical depths concerning individual features refer to the AGN component only, and are evaluated considering the 
observed (i.e. not corrected for reddening) AGN continuum from the best fits. 
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Moreover, according to a recent spectroscopic survey of young stellar objects, both water 
ice and HAC seem to play a deficient role among the carriers of the most enigmatic 5-8 fim 
depressions (Boogert et al. 2008). This could provide a physical explanation also to the 
spectral peculiarities observed in the three anomalous cases, in which the deviations from 
the standard ~6 /xm and 6.85 /xm profiles give rise to a formless band that undermines our 
approach. 

8.4 Absorption versus reddening 

An accurate chemical analysis goes beyond the aims of the present work, anyway a further 
line of investigation can be pursued, involving the study of possible correlations between 
the optical depth of the most prominent absorption features and the reddening of the dust 
continuum. In Figl9] some representative ULIRG spectra have been arranged in a slope se- 
quence with decreasing steepness with respect to the observed AGN component. Different 
colours have been used in order to flag the absorption degree and avoid ambiguities in the 
following considerations. Through a simple comparison among the adjacent spectra, it is 
clear that the intensity of the absorption features varies considerably within the same range 
of continuum slope. The effects of extreme absorption are evident e.g. in the black spectra 
(IRAS 20100-4156, IRAS 00188-0856 and UGC 5101), all of which also display strong 
aromatic emission longward of 6 fim. At shorter wavelengths the observed flux density is 
unduly intense to share a common origin with the PAH features: this can be easily inferred 
by considering the shape of our SB template. If such emission is to be interpreted as the 
AGN-related hot dust continuum, one is in the presence of really contrasting signatures, 
since the aromatic nature of both the 6.2 and 7.7 peak is beyond dispute. Conversely, 
in cases like the blue spectrum (IRAS 01298—0744) on the top of FigJHl the possibility that 
the tentative 6.2 fim feature is actually a ridge in the absorbed continuum mimicking the 
PAH emission can not be ruled out. Moreover the supposed 7.7 /im companion appears to be 
broad and roundish, as usually observed in the AGN-dominated spectra (as also the green 
ones, i.e. IRAS 16156+0146 and IRAS 08572+3915): this is rather a narrow window of unab- 
sorbed continuum, that immediately declines into the silicate trough. A possible explanation 
to the anomalous spectra can be the coexistence of a buried AGN component of moderate 
luminosity with a powerful starburst. In this scenario the direct AGN emission is quenched 
by a dense screen, producing deep absorption bands whose carriers are unknown (perhaps 
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Figure 9. From top top bottom, the 5-8 /xm spectra of the following sources are plotted in a slope sequence with decreasing 
steepness: (1) IRAS 01298-0744, (2) IRAS 20100-4156, (3) IRAS 16156+0146, (i) IRAS 00188-0856, (5) IRAS 08572+3915, 
(6) IRAS 21219-1757, (7) IRAS 12127-1412 and (8) UGC 5101. Colours are used to denote the consequence of absorption 
{red, green, blue and black in order of growing significance). This contrast shows the wide variety in both reddening and 
absorption characteristic of ULIRGs. 

various ices and HAC, see above). The sharp disappearance of the hot dust continuum is 
justified if the blending reaches up to the 9.7 yum sihcate feature. A hint in this direction 
is provided by the comparison between the contiguous spectra of IRAS 21219—1757 {red, 
unabsorbed) and IRAS 12127—1412 {blue): both these sources can be considered as pure 
AGN at 5-8 /im. In the former we observe a featureless and uniform continuum, which is 
a reasonable candidate to outline the true hot dust emission of the latter source as well. In 
this case, however, the measured 8 /im fiux density is well below that expected from the 
extrapolation of the putative continuum slope. Hence it possibly represents only a saddle 
point between the broad ~7 and 9.7 /im absorption bands. 

The qualitative considerations above are of course insufficient to claim that very unusual 
physical conditions are responsible for the emergence of anomalous spectral shapes and 
absorption/reddening polarities. In order to properly investigate these aspects and charac- 
terize the anomalous sources, a detailed analysis is needed from if-band up to ~20 /^m. This 
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would make it possible to address the problems connected to the uncertain location of the 
continuum. For a few of our targets this spectral range is nearly covered in its entirety, and 
interesting indications come from L- and M-band studies in which strong aliphatic and CO 
absorption features are detected: the red wing of the 4.65 /xm CO profile is indeed entering 
the very beginning of the Spitzer-IKS window in several absorbed objects of our sample. 
None the less a joint modelling of the 2-20 fim SED is at present frustrated by the disparity 
in signal to noise and the scatter in flux calibration among the different ground-based and 
space facilities involved in the observations. 

Summarizing, there seems to be no clear correlation between the degree of continuum red- 
dening and the optical depth of individual absorption features. From a quantitative point of 
view this has already been established by Sani et al. (2008) in an extensive L- and M-band 
study of five bright ULIRGs. The presence of anomalous sources in our sample challenges 
even the possibility of a step-wise correlation, which is instead found shortward of 5 fim: 
neither the detection of broad absorption features implies invariably a sizable reddening, nor 
a flat continuum is regularly unabsorbed. This should not be surprising, since ULIRGs have 
a very disturbed morphology and large variations in chemical composition and geometrical 
structure of the dust are naturally expected. Yet it is possible to reconcile the observation of 
significant absorption features alongside apparent lack of reddening without invoking special 
dust distributions and covering factors. In Fig {10] the optical depths of the main absorption 
features are plotted against the optical depth to the hot dust continuum. All these values 
are also included in Tabj3l By focusing on the confident entries only, the positive correlation 
between absorption and reddening is tentatively recovered. This suggests that modifications 
of the extinction law are possibly involved. 

8.5 The extinction law 

Conforming to our interpretation, both reddening and deep absorption are related to the 
presence of a compact power source, and must be due almost exclusively to cold material in 
the surroundings of the AGN component. We therefore argue that a reddened continuum is 
not observed in the anomalous sources because the actual extinction law applicable to their 
AGN component has a softer wavelength dependence. Within the power-law prescription 
Tx oc A^^ ''^ that has been adopted after Draine (1989), the ratio between the values of the 
optical depth at the red and blue ends of our fitting region is Tred/Tuue —0.44; this allows 
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Figure 10. Opacity of absorption features versus continuum opacity, inferred from the reddening of the AGN component. 
Prom top to bottom, the three panels show the qualitative correlation between the optical depth of the main observed features 
(ascribed to water ice at ~6.0 /im and hydrogenated amorphous carbons at ~6.85 /im and ~7.25 fim) and the optical depth to 
the hot dust continuum. In each plot, the blue crosses refer to well-established measurements of both quantities, while the other 
points render in a different colour code all the cases involving upper limits. Only the 50 AGN detections are considered, since 
no absorption is found in pure starbursts. In a few cases the observed AGN component is too faint and no sort of indication can 
be evinced about additional features: for this reason, the number of effective entries is 48, 43 and 43, respectively. The values 
of the optical depth for the three anomalous sources are computed assuming a flat AGN intrinsic continuum and a quasi-grey 
extinction law (see the discussion in the text) and are provided as upper limits on the vertical axis, since the identification 
of the single features is unclear. If we take into account the secure entries only, we evince that no strict proportionality holds 
between reddening and absorption, nevertheless they positively correlate. This evidence is not controverted by the other points, 
whose constraints are very limited. We conclude that the appearance of broad absorption coupled with reddening deficit is 
confidentially due to a flat extinction law. 



US to reproduce the steepest continuum gradients. Yet, according to the observations, the 
extinction curve can be safely described by a power law only up to ~4-5 fim, and becomes 
rather controversial in the 5-8 fim range. For a long time, in absence of an adequate coverage 
of this spectral region with scientific data, the near-IR power-law trend was supposed to hold 
up to the 9.7 fim silicate bump. Only in the last decade, in the wake of the observations 
towards the Galactic Centre (Lutz et al. 1996), evidence has grown that the extinction curve 
can also undergo a flattening longward of ~4 fim, without giving rise to a clear minimum 
before the silicate feature (e.g. Indebetouw et al. 2005; Nishiyama et al. 2009, and references 
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therein). Since we are applying extinction to a compact nuclear source, this flat variant rep- 
resents another reasonable choice. Moreover all the anomalous sources (regardless of their 
absorption features) have a very steep L-band continuum, suggesting that the AGN com- 
ponent is actually reddened and that some change in the extinction pattern may occur at 
longer wavelengths. 

In order to test this scenario, we have performed our AGN/SB decomposition once again, 
this time assuming the analytical form of Chiar & Tielens (2006) for the wavelength depen- 
dence of the optical depth. The intrinsic slope of the AGN template remains ragn=^-5- The 
alternate extinction law is really quasi-grey in our working range, since Tred/Tuue —0.90 (and 
'Tred/'TQ —0.96). This eutails a degeneracy between the AGN flux amplitude and its optical 
depth, in spite of which we are able to reach important conclusions. The overall results 
can be summarized in terms of subclasses, (a) A fiat extinction can be safely ruled out 
for eight sources, whose observed continuum is too steep to be reproduced within the new 
framework. This group encompasses all the objects with aboi > 0.25 and rg > 1 according 
to the previous estimate; the only exception is IRAS 15250+3609. (b) The objects whose 
classification is unchanged reach up to 30, by combining 19 out of 21 pure starbursts (save 
IRAS 02411+0353 and IRAS 16487+5447, that may contain a faint AGN component) and 
11 out of 14 composite sources hosting a non-extinguished AGN. (c) The remaining three, 
along with IRAS 07598+6508, IRAS 19254-7245, the above-mentioned IRAS 16487+5447 
and other four mildly reddened sources, have a non-degenerate fitting output, (d) Finally, 
there are 23 objects (including the anomalous ones) that can be fitted with the modified 
extinction but are affected by the normalization/absorption degeneracy. It has been possible 
to remove such degeneracy in 20 cases, by pegging the value of the optical depth after our 
previous L-band studies. 

The consequence of a quasi-grey extinction have been tested by computing the new values of 
R and and checking their re- arrangement with respect to the former best fit of the R-a^ 
relation. The outcome is shown in FigHHand implies two interesting results: 1) a fiat extinc- 
tion law seems suitable to describe the dusty environment surrounding the AGN component 
in many ULIRGs; 2) in a statistical sense, the assumption of a quasi-grey extinction does 
not alter our global results, once the optical depth is reliably estimated through additional 
constraints. In conclusion, the large size and the spectral variety of our sample suggest that 
the AGN component inside a ULIRG can not be described in terms of a universal 5-8 /xm 
extinction pattern: both a power-law and a soft wavelength dependence of the optical depth 
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Figure 11. R-ag plot obtained by assuming a quasi-grey extinction law for the AGN component. The new positions of the 
sources are compared with those obtained with a power-law extinction: the orange solid line is indeed the best fit of the R-a^ 
relation from Eq. Different colours and symbols are used according to the following: black crosses for sources confirmed to 
host a non-extinguished AGN component, if any; blue rings for sources with a definite fit in the new scenario; red diamonds 
for sources whose normalization has been fixed through their L-band spectra; green triangles for sources in which the removal 
of the normalization/absorption degeneracy is tentative, due to the lower quality of the L-band spectra; magenta points for 
sources appearing as pure starbursts in the L-band. 

appear to be involved. In addition to a different dust amount and composition along the 
line of sight, this is possibly due also to geometrical effects, such as the orientation of the 
obscuring torus itself and/or the dumpiness of the circumnuclear star- forming cores. It is 
however important to keep in mind the limitations of our approach: our description does 
not take into account the possible effects of radiative transfer, that can be essential even at 
5-8 yum. In particular, the deviations from a power-law extinction could be related to a lesser 
accuracy of the screen approximation, which is instead correct for the Galactic Centre. This 
issue represents a critical challenge for any future study of the dust properties and near- and 
mid-IR extinction in AGN and ULIRGs. 

8.6 Trend with bolometric luminosity 

Once the relative AGN/SB fraction of ULIRG luminosity has been estimated, one of the main 
issues to investigate is the possibility of a larger AGN contribution at higher luminosities. 
Interestingly, no clear trend in the detection of [Ne v] A 14. 32 as a function of the total 
IR luminosity has been discovered by Farrah et al. (2007). This is in contrast with several 
studies based on optical spectroscopy, which instead suggest the existence of a growing 
trend. This claim finds the first convincing confirmation with the histograms of Veilleux 
et al. (1999a), indicating that AGN-like systems (i.e. those classified as Seyferts, without 
considering the possible integration of LINERs) represent only a small percentage among 
IR galaxies below Ljr < IO^^Lq, but increase to become the most populated subset in the 
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luminosity range IO^^'^Lq < Ljji < IO^^'^Lq. Moreover virtually all the hyper luminous IR 
sources (HLIRGs, Ljr > IO^^Lq) seem to harbour a powerful AGN, at least up to 2; < 1 
(Rowan- Robinson 2000). The key point is that such a result still has to be read in terms of 
AGN detectability, and not necessarily of larger AGN contribution. In a skeptical perspective 
even this detectability trend could be contrived by selection effects: as we have already 
pointed out, a considerable amount of sources that are optically classified as H ll regions 
actually harbour obscured but non-negligible AGN components. Furthermore, a general 
consensus has grown in the last years about the sharp decrease in the fraction of absorbed 
AGN with X-ray luminosity (Ueda et al. 2003; Hasinger 2008). At mid-IR wavelengths a 
clear AGN preponderance at high luminosity is found by Tran et al. (2001), with starbursts 
and AGN dominating the energy output respectively at the lower and higher ends of the 
ULIRG range. It is then worthwhile checking how our results fit into this question, even 
though our coverage of bolometric luminosity is limited. 

Assuming the standard divide of log^Lju/ Lq) = 12.3, the population of the two subsets 
consisting of the sources whose luminosity is below and above this threshold is 55 and 16, 
respectively. We recall that simply by solving the ambiguity concerning LINERs we are able 
to detect a much larger number of AGN (50 out of 71 ULIRGs) with respect to optical 
studies. The positive detections involve ~2/3 of the faint sources and ~3/4 of the bright 
ones. Such fractions are too similar to hint at a clear increase. Nevertheless by computing 
the AGN contribution aagn to both bins of luminosity we obtain two well-separated values: 

_ 18.2tli (log(L,fi/LQ) < 12.3) 

an 



'■agn 



30.9tl:^ (log(Lffi/Lo) > 12.3) 
where a and L stand for the previous aboi and Lj^. This is robust evidence in favour of 
the growing trend with bolometric luminosity of the AGN content inside ULIRGqj. An 
illustration of this effect, that will be addressed in detail in a forthcoming paper (Nardini et 
al., in preparation) by enlarging the statistics at the highest luminosities (see also Imanishi 
2009), is provided in FigUl 

* There is indeed a difference in the cosmological parameters between this work and that of Kim & Sanders (1998), on which 
the mentioned investigations of the 1 Jy sample at optical wavelengths are based. This slightly affects the estimate of the 
bolometric luminosities, with an upward shift of ~0.1 dcx at the highest luminosity distances: about ten sources enter the 
higher luminosity bin. We have taken into account this effect, and found that the trend is both qualitatively and quantitatively 
confirmed. 
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Figure 12. Relative number counts witiiin tiio different bins of AGN contribution. The two distributions correspond to the 
populations of sources with \og{Ljji/ Lq) < 12.3 {cyan) and \og{Ljji/ Lq) > 12.3 {shaded blue). 

In conclusion, considering only the 63 sources from the 1 Jy sample, AGN are responsible 
for 23.llg g per cent of the ULIRG luminosity in the local Universe. We underline that this 
estimate is not biased due to selection effects. Concerning the seven missing sources which 
undermine the sample completeness, by simply assuming the average properties of their 
optical class for the five type 2 Seyferts (i.e. ai^oi —0.29) and a SB-dominated nature for 
the two unclassified objects we are able to confirm the value provided above. Only in the 
extreme case of null AGN significance for all the seven sources the global contribution of 
accretion processes to IR activity should be reduced by ~2 per cent. The same amplitude in 
the opposite direction is predicted after correcting for a possible bias towards star formation, 
connected to the cold 60 /im selection. Not even the inclusion of the eight southern ULIRGs 
in the computation affects our results, since they form a complete and fully representative 
subset. 

The recent and comprehensive work of Veilleux et al. (2009) investigates with great detail the 
connection between local ULIRGs and PG quasars; the contribution of nuclear activity to the 
bolometric luminosity of ULIRGs is obtained by means of six independent methods based on 
the Spitzer-IR^ spectra. Fine-structure line ratios, aromatic features and continuum colours 
have been employed: a comparison among the six different estimates provided for the single 
sources gives a good idea of the uncertainties involved, that sometimes can be pretty large. 
The ensemble average hints at a slightly higher AGN contribution (~35-40 per cent) with 
respect to the findings of our study and the previous literature. This could be due to a 
possible bias in Veilleux et al. (2009), related to the choice of PG quasars as zero points 
for the nuclear activity within ULIRGs. However, by simply restricting to the 48 common 
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sources, our estimate increases to ~30 per cent as well, and the results are definitely in good 
agreement. 

9 CONCLUSIONS 

In this paper we have presented a 5-8 /xm spectral study of a large and representative 
sample of local ULIRGs, based on Spitzer-UXS) observations. We have explored the role of 
black hole accretion and vigorous star formation as the power source of the extreme IR 
activity, assessing the gravitational and stellar origin of ULIRG emission by modelling their 
5-8 //m spectra through a couple of AGN and SB templates. In this wavelength range, in 
fact, the spectral properties of active galactic nuclei and star bursts are widely different, and 
only moderate dispersion is found within each class. It is then possible to accomplish a sharp 
characterization of both the AGN and SB components, and disentangle their contribution 
to the observed emission. Our method has proven to be successful in unveiling intrinsically 
faint or heavily obscured AGN components, leading to a total of 50 convincing detections 
out of the 71 sources in our sample. Such a large AGN detection rate can be achieved only 
by resorting to the most effective multiwavelength diagnostics. As a consequence, AGN turn 
out to be very common among ULIRGs in the local Universe, even if they are usually not 
significant as contributors to the global energetics of the host galaxy when compared to the 
starburst. 

A simple analytical model, in fact, allows us also to obtain a rehable estimate of the AGN/SB 
contribution to the overall energy output of each source. For the sake of a gross classification 
we can define the following one third rule: assuming a significance threshold of 1/3 in the 
AGN/SB bolometric contribution, i.e. a^oi — 0.25, we find that one third of ULIRGs harbour 
a sizable AGN, one third harbour a negligible AGN, and one third harbour no AGN at all. 
Star formation is then confirmed as the dominant energy source underlying the ULIRG 
activity. In more quantitative terms, AGN account for ~23 per cent of the global ULIRG 
emission, approaching and possibly exceeding ~30 per cent at higher luminosities. This 
increasing trend is clearly brought out: the average AGN weight rises from 18 to 31 per cent 
above the luminosity threshold of Lm — lO^^-^L©. 

Thanks to the large size of our sample we have also been able to test the properties of the 
obscuring material, in particular of the extinction law applicable to the AGN-related hot 
dust continuum. Observations reveal a lack of correlation between the continuum reddening 
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and the presence of deep absorption features, suggesting that the extinction of the AGN 
component in a ULIRG environment is not universaL Both a power-law and a quasi-grey 
behaviour of the optical depth as a function of wavelength are necessary to account for 
the emission of different objects and seem to be involved among ULIRGs. Large differences 
in the chemical composition and geometrical structure of the dust are supported by the 
spread in the intensity (and sometimes also the shape) of the most frequent absorption 
profiles. Although our method is not suitable to fully investigate these aspects, it is extremely 
powerful for global classification and can be applied whenever spectra of sufficient quality 
are available. Moreover it can be turned into a photometric method to study fainter sources, 
using diagnostics such as the bolometric correction and the continuum slope, which can be 
measured efficiently with the upcoming IR facilities. 
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APPENDIX A: ADDITIONAL NOTES AND SPECTRAL FITS 

With the exception of IRAS 08559+1053, that was observed in the early stage of the mission 
(before the start of the regular scientific operations) and is analyzed for the first time in 
the present work, the low- and/or high-resolution 5'pzfeer-IRS spectra of the sources in our 
sample have already been shown and widely discussed by Armus et al. (2007), Farrah ct 
al. (2007), Imanishi et al. (2007), and Vcillcux ct al. (2009). In the following we add some 
additional comments about a few interesting objects, and show the full results of our spectral 
decomposition. 



Al The nature of ARP 220 

In spite of being the nearest (at ~80 Mpc) and best-studied ULIRG, ARP 220 is still 
puzzling and controversial with respect to its energy source. The multi wavelength SED of 
this interacting system has been often employed as a template to constrain the emission of 
star-forming galaxies at high redshift, but the pure SB scenario is actually challenged by 
many independent pieces of observational evidence, and the alleged action of an extremely 
obscured and elusive black hole in the western nucleus is widely debated. Such a possibihty 
is strengthened by recent interferometric observations, outlining a compact and massive ring 
of hot dust that straightly recalls the toroidal structures around AGN (Downes & Eckart 
2007, and references therein for additional hints to the presence of an AGN inside ARP 220). 
It is somewhat arbitrary to assess the intrinsic luminosity of this core; the western nucleus 
is anyway the brightest one by a factor ~3 (Soifer et al. 1999) and the AGN component, 
if confirmed, would be absolutely relevant. Our diagnostic method discloses a contribution 
to the bolometric luminosity of the whole system arising from a highly-obscured AGN, that 
reaches up to 20 per cent. This value should be read with some caution, also due to the 
peculiar properties of the aromatic emission in ARP 220 (e.g. the suppression of the 6.2 //m 
feature with respect to the 7.7 //m one) that may alter the continuum decomposition based 
on AGN/SB templates; yet it represents a stirring indication of magnitude. The matter 
is fully open, and further incentive is provided by Sakamoto et al. (2008), whose precise 
discussion about the nature of the compact luminous source located in the western nucleus, 
both in terms of an AGN and of a collection of overlapping super star clusters, makes the 
latter explanation appear as slightly strained. 
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Table Al. Flux and equivalent width (with respect to the AGN component) of the 7.65 fim [Ne vi] fine-structure line. The 
uncertainties on both quantities range from ~5 to 30 per cent. 



Source 


[Ne vi] Flux 


EW 


(IRAS) 


(ergs s"-*- cm~^) 


(/im) 


04103-2838 


1.0 X 10^13 


0.024 


05189-2524 


1.1 X 10^13 


0.008 


12072-0444 


2.3 X 10"" 


0.011 


13454-2956 


1.5 X 10"" 


0.040 


17179+5444 


3.8 X 10-" 


0.057 



A2 The iron Ka line in IRAS 04103-2838 

IRAS 04103—2838 is optically classified as a LINER and is a warm ULIRG according to 
its IRAS colours (/25//60 = 0.30). We are able to unveil an AGN component at 5-8 fim, 
whose contribution to the bolo metric luminosity is estimated as ~6 per cent. The presence 
of an AGN is confirmed by the detection of a broad Ka line at ~6.4 keV, during a recent 
XMM-Newton observation (Teng et al. 2008). Due to the large equivalent width (~1.65 keV) 
of this iron feature the AGN is supposed to be Compton-thick, anyway it does not seem to 
dominate the total energy output of the source. 



A3 The 7.65 fim fine-structure line 

IRAS 17179+5444 exhibits a prominent emission at ~7.65 fim, blended with the adjacent 
PAH feature, that we interpret as the unresolved [Ne vi] fine-structure line. Due to its 
ionization energy of 126 eV this feature is connected to the hardest nuclear activity, and in- 
deed it is convincingly observed also in other three Seyfert galaxies (IRAS 05189—2524, 
IRAS 12072-0444 and IRAS 13454-2956) plus IRAS 04103-2838. Fluxes and equiva- 
lent widths are listed in Tab lAlt and are in good agreement with the values provided by 
Veilleux et al. (2009) for their positive detections. We note that our identification of this 
line also in IRAS 04103—2838 represents the first detection involving a non-Seyfert ULIRG. 
In IRAS 17179+5444 the strength of [Ne Vl] is anomalously high, as evinced through both 
a visual inspection and a comparison with the [Ne v] lines at longer wavelengths (Farrah et 
al. 2007). This is a very interesting case, since after a separate check of the single spectra 
from each nod position the line is confirmed to be genuine. 
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A4 The case study of NGC 6240 

As mentioned, many X-ray observations suggest the possible presence in NGC 6240 of a pair 
of supermassive black holes hidden behind large column densities of screening material. A 
detailed IR study pivoted on the ^piteer-IRS spectrum and integrated with photometric data 
from near- to far-IR is presented by Armus et al. (2006). As a result of the SED modelling 
a hot dust component at ~700 K is found, accounting for ~3.5 per cent of the bolometric 
luminosity and pointing to the presence of a buried AGN. The AGN contribution is again 
estimated as ~3-5 per cent through high-ionization lines, but can raise up to ~20-24 per 
cent if correcting the [Ne v]/[Ne ll] ratio for the extinction derived from the hard X-ray data. 
Our estimate of ahoi (4.5-11.5 per cent) is then fully reliable, in spite of the great simplicity 
of our method as opposed to a multicomponent fit of the whole IR SED, which is anyway 
necessary to fully characterize an individual source. 

A5 Spectral fits 

The following panels illustrate the results of our AGN/SB decomposition at 5-8 jim.. Along 
with the data points {green filled circles) and the best fit model [black solid line), both 
contributions are plotted: the SB component [red dot-dashed line) and the observed AGN 
continuum [blue long-dashed line). 
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